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Abstract 

Two-dimensional  regular  and  Mach  reflections  have  been  studied  in  the  Mach 
4.96  dual-solution  domain  for  a  25°  and  26°  double-fin  inlet.  The  steady-state  com¬ 
putational  Mach  and  regular  reflections  were  subjected  to  magnetogasdynamic  forces 
to  determine  whether  these  forces  could  be  used  as  a  possible  flow  control  mechanism. 
The  numerical  code  employed  for  this  research  solved  the  inviscid  Euler  equations 
with  added  source  terms  for  the  ponderomotive  force  and  accompanying  energy  in¬ 
teractions.  The  25°  regular  reflection  was  determined  to  be  extremely  sensitive  to  a 
decelerating  Lorentz  force.  Transient  application  of  the  force  led  to  the  transition  of 
the  regular  reflection  to  a  Mach  reflection,  increasing  the  total  pressure  losses  and 
decreasing  the  compression  ratio.  Sustained  application  of  the  force  resulted  in  inlet 
unstart.  An  accelerating  Lorentz  force  was  also  examined  with  the  goal  of  transition¬ 
ing  the  26°  Mach  reflection  to  a  more  efficient  regular  reflection.  The  location  of  the 
accelerating  force  and  the  parameters  governing  its  magnitude  were  examined.  Such 
forces  push  the  Mach  reflection  back  to  a  more  stable  location  and  reduce  the  Mach 
stem  height.  For  the  interaction  parameters  considered,  fully  regular  reflections  were 
not  obtained.  However,  the  accelerating  Lorentz  force  proved  capable  of  increasing 
the  total  pressure  recovery  and  the  static  pressure  compression  beyond  the  regular 
reflection  values. 
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MAGNETOGASDYNAMIC  FLOW  CONTROL 


OF  A  MACH  REFLECTION 


1.  Introduction 

Starting  with  the  first  flight  of  a  powered  aircraft  on  17  December  1903  by 
Orville  and  Wilbur  Wright,  mankind’s  continuing  desire  to  fly  ’’faster  and  higher”  [1] 
has  pushed  the  limit  of  aircraft  and  spacecraft  design.  These  aspirations  have  led  to 
many  new  performance  achievements,  from  the  first  hypersonic  flight  of  the  modified 
two-stage  V-2  rocket  in  1949  [1]  to  the  Space  Shuttle  Columbia’s  first  manned/winged 
vehicle  re-entry  flight  [12] .  Figure  1.1  depicts  a  few  of  the  milestones  of  high-speed 
flight  in  the  past  century.  These  milestones  have  helped  pave  the  way  in  hypersonic 
flight,  leading  towards  the  ultimate  goal  of  sustained  hypersonic  flight  whereby  one 
could  fly  from  New  York  to  Tokyo  in  as  little  as  two  hours[l].  The  cost  of  putting  a 
payload  in  space  could  be  reduced  by  l-to-2  orders  of  magnitude  with  a  single-stage- 
to-orbit  (SSTO)  hypersonic  aerospace  vehicle. 

Hypersonic  flight  can  be  achieved  with  rocket  or  air-breathing  propulsion. 
Rocket  propulsion,  both  solid-  and  liquid-fuelled,  is  the  current  means  for  hypersonic 
travel.  Over  the  past  fifty  years,  many  missiles,  aircraft,  and  spacecraft  have  been 
propelled  to  hypersonic  velocities,  greater  than  Mach  5,  by  rocket  engines.  However, 
the  objective  of  sustained  hypersonic  flight,  manned  or  unmanned,  drives  propulsion 
designers  toward  an  air-breathing  propulsion  source.  Air-breathing  engines  use  at¬ 
mospheric  oxygen  instead  of  carrying  an  oxidizer.  Subsequently,  the  payload  and/or 
the  range  of  the  hypersonic  vehicle  increases.  Therefore,  hypersonic  transports  could 
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Figure  1.1  Time-line  of  major  hypersonic  and  space  flight  events [1,  12,  16,  21,  24] 


carry  a  greater  payload  for  longer  distances,  and  hypersonic  missiles  could  project 
larger  warheads  over  a  greater  range  when  the  internal  oxidizer  is  removed. 

One  of  the  air-breathing  propulsion  types  studied  for  hypersonic  flight  is  the 
supersonic  combustion  ramjet  (SCRAMjet).  Figure  1.2  depicts  the  National  Air  and 
Space  Administration’s  (NASA)  X-43A[21],  The  X-43A  is  an  unmanned  hypersonic 
vehicle  which  uses  a  SCRAMjet  to  fly  between  Mach  =  7  and  10.  The  vehicle  is 
boosted  to  Mach  7  by  a  Pegasus  solid-rocket  motor  and  then  separates  from  the 
rocket  and  ignites  the  SCRAMjet.  The  first  X-43A  flight  was  lost  in  June  of  2002 
due  to  problems  caused  by  the  solid-rocket  booster. 

The  X-43  demonstrates  many  of  the  key  features  of  a  SCRAMjet  engine.  The 
initial  compression  and  slowing  of  the  flow  is  due  to  the  bottom  surface  of  the 
aircraft  as  depicted  in  Fig.  1.2.  This  surface  acts  as  a  compression  ramp  creating 
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Underside  of  Aircraft 


Variable  Angle 
Compression  Ramps 

Figure  1.2  Schematic  of  the  X-43A  hypersonic  vehicle  [NASA] 


an  oblique  shock  ending  on  the  lower  lip  of  the  inlet  [2].  The  flow  then  goes  through 
a  complex  set  of  oblique  shocks  created  by  the  variable  angle  compression  ramps  in 
the  inlet  of  the  engine.  These  ramps  are  adjusted  to  the  angle  required  to  create  the 
proper  compression  for  the  flight  condition.  The  supersonic  mixing  and  combustion 
of  the  fuel  occurs  in  the  combustor  section.  The  resulting  high-temperature  gases 
are  expanded  through  the  nozzle  and  ejected  out  the  rear  of  the  engine [10].  The 
exhaust  gases  are  further  expanded  to  freestream  pressure  using  both  the  bottom 
surface  of  the  aircraft  behind  the  engine  and  the  pressure  difference  of  the  freestream 
flow [2].  This  sequence  of  events  is  controlled  completely  by  the  variable  ramps  and 
the  injection  of  fuel. 

The  complex  shock  wave  interactions  present  in  internal  flowfields  influence 
the  design  and  performance  of  SCRAMjet  engines.  The  shock  waves  interact  with 
each  other  and  the  turbulent  boundary  layers  on  the  adjacent  walls  creating  com¬ 
plex  shock-shock  and  shock-boundary  layer  structures.  This  highly  non-linear  phe¬ 
nomenon  can  induce  separation  in  the  boundary  layer  and  extreme  deformation  in 
the  flowfield.  These  effects  lead  to  total  pressure  losses,  reduced  performance,  and 
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possible  engine  unstart  where  the  shock  system  is  expelled  from  the  inlet.  These 
potential  problems  drive  researchers  and  designers  to  attempt  to  understand  and 
control  the  complex  interaction  phenomena. 

The  complicated  geometry  of  a  SCRAMjet  inlet  has  been  modelled  recently  as 
shown  in  Fig.  1 .3 [IT,  20,  29,  30].  This  model  is  composed  of  a  flat  plate,  represent¬ 
ing  the  underside  of  the  aircraft  compression  ramp  in  Fig.  1.2,  and  perpendicular 
shock  generators  on  each  side  of  the  flat  plate.  The  oblique  shock  generated  by  the 
underside  of  the  aircraft  is  taken  into  account  by  the  freestream  conditions  entering 
the  inlet,  while  the  flat  plate  simply  grows  a  boundary  layer  which  induces  shock¬ 
boundary  layer  interactions.  The  shock  generators  can  be  adjusted  to  any  angle  6  to 
change  the  angles  of  the  oblique  shocks,  /?,  and  the  corresponding  flow  compression 
of  the  inlet. 

1.1  Regular  and  Mach  Reflections 

When  Moo  =  4.96,  there  exists  an  essentially  two-dimensional  (2-D)  region  for 
this  double-fin  configuration  beginning  27  boundary  layer  thickness  above  the  flat 
plate [29].  In  order  to  investigate  this  region,  a  2-D  domain  is  created  from  a  single 
slice  of  the  three-dimensional  (3-D)  model  in  Fig.  1.3.  The  slice,  shown  in  Fig.  1.4, 
is  used  to  calculate  the  interactions  of  the  oblique  shocks  without  the  influence  of 
the  boundary  layer  of  the  flat  plate.  In  agreement  with  theory,  experimental  and 
computational  results  show  that  either  a  regular  reflection  or  a  Mach  reflection  can 
exist  at  these  conditions  when  6  is  between  22°  and  28°  [13,  14,  17,  20,  29,  30]. 

A  regular  reflection  (RR)  is  characterized  by  oblique  shocks  off  the  leading 

edges  of  the  fins  at  the  angle  specified  by  oblique  shock  theory.  The  shock  angle,  /?, 

is  a  function  of  the  freestream  Mach  number  and  the  compression  angle,  9.  Figure  1.5 

shows  the  RR  generated  by  a  double  fin  configuration  with  9  =  25°  and  M^  =  4.96. 

At  the  intersection  of  the  primary  oblique  shocks  generated  by  each  fin,  the  shocks 

reflect  in  a  regular  fashion.  Expansion  fans  form  at  the  throat  of  the  inlet  and 

interact  with  the  reflected  oblique  shocks  inducing  curvature  towards  the  walls.  The 
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Figure  1.3 


Figure  1.4 


Schematic  of  the  double-fin  generated  crossing  shock-wave  turbulent 
boundary-layer  interaction 


Schematic  of  the  2-D  slice  from  the  double-fin  configuration  turbulent 
boundary-layer  interaction 
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oblique  shocks  continue  to  reflect  off  the  walls  and  each  other  as  the  flow  continues 
through  the  inlet  until  they  exit  the  channel  as  illustrated  by  Fig.  1.5. 


Figure  1.5  Pressure  contours  from  a  regular  reflection  for  =  4.96,  2-D  solu¬ 
tions,  25°  x  25° 

The  Mach  reflection  (MR)  for  the  equivalent  conditions  has  oblique  shocks 
coming  off  the  leading  edges  of  the  fins  identical  to  the  regular  reflection.  However, 
a  normal  shock  (Fig.  1.6)  is  formed  between  the  two  oblique  shocks  instead  of  the 
oblique  shock  to  shock  reflection.  The  size  of  the  normal  shock  is  a  is  determined 
by  the  freestream  Mach  number  and  fin  angles.  This  normal  shock  is  also  known  as 
the  Mach  stem. 

The  region  behind  the  MR  is  characterized  by  a  slip- line  between  the  subsonic 

flow  behind  the  normal  shock  and  the  supersonic  flow  behind  the  reflected  oblique 

shocks  as  shown  in  Fig.  1.7.  The  large  region  of  subsonic  flow  decreases  the  inlet 
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Figure  1.6  Pressure  contours  from  a  Mach  reflection  for  =  4.96,  2-D  solutions, 
25°  x  25° 

efficiency  and  increases  the  total  pressure  losses.  This  Mach  reflection  could  occur 
in  the  inlet  due  to  changes  in  the  performance  of  the  propulsion  system  at  off-design 
flight  conditions. 

Whether  the  flow  develops  a  RR  or  a  MR  is  dependent  on  the  pressure  rise 
across  the  oblique  shocks[13,  14].  The  von  Neumann  angle,  9n,  is  the  wedge  angle 
for  a  given  freestream  condition  where  the  pressure  rise  across  the  initial  and  re¬ 
flected  oblique  shocks  is  equal  to  the  pressure  rise  across  a  single  normal  shock.  The 
detachment  angle,  6d,  is  the  wedge  angle  for  a  given  freestream  condition  when  the 
initial  oblique  shock  waves  are  too  strong  for  the  reflected  oblique  shock  waves  to 
remain  attached  to  each  other.  Ramp  angles  less  than  9 ^  form  regular  reflections 
and  ramp  angles  greater  than  9 d  form  Mach  reflections.  The  range  of  wedge  angles 
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Figure  1.7  Isobars  of  Mach  number  from  a  Mach  reflection  for  =  4.96,  2-D 
solutions,  25°  x  25° 

between  6n  and  0(i  for  a  given  freestream  Mach  number  is  known  as  the  dual-solution 
domain  where  both  a  RR  and  a  MR  can  exist [13,  14]. 

The  inlet  ramp  angles  are  adjusted  to  a  specific  6  to  produce  a  set  of  oblique 
shocks  with  the  proper  inlet  compression  ratio.  The  required  compression  ratio  is 
based  on  the  flight  condition  and  the  throttle  position.  A  MR  produces  a  lower  com¬ 
pression  ratio  resulting  in  reduced  thrust.  In  addition,  normal  shocks  have  greater 
losses  in  total  pressure  across  them  than  oblique  shocks.  Therefore,  the  MR  condi¬ 
tion  has  higher  total  pressure  losses  and  lower  inlet  efficiency  than  the  corresponding 
RR  as  shown  in  Fig.  1.8.  These  large  total  pressure  losses  associated  with  the  nor- 
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mal  shock  in  the  Mach  reflection  can  lead  to  an  unstart  of  the  engine  inlet.  [9]  The 
consequences  of  engine  unstart  include  a  loss  in  the  total  pressure  recovery,  dramatic 
increases  in  drag,  and  potential  loss  of  aircraft  control. 


P, 

617.887 

477.134 

336.382 

195.63 

54.8773 


Figure  1.8  Mach  and  Regular  reflection  total  pressure  comparison  for  =  4.96, 
2-D  solutions,  25°  x  25° 

1.2  Previous  Research 

Extensive  research  has  been  performed  by  several  researchers  on  this  dual¬ 
solution  domain[17,  20,  29,  30,  34],  The  research  includes  numerical  simulations 
and  experiments  for  both  two-  and  3-D  configurations.  The  3-D  computations  and 
experiments  explore  the  shock-shock  and  shock-boundary  layer  interactions,  and  how 

these  interactions  affect  each  other.  Although  the  3-D  research  has  been  conducted 
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in  the  dual-solution  domain,  these  simulations  and  experiments  have  not  been  able 
to  demonstrate  a  regular  reflection,  resolving  a  Mach  reflection  in  every  case  [17, 
18,  20,  29,  30].  This  is  attributed  to  the  boundary  layer’s  effect  on  the  oblique 
shock  interaction  and  the  accompanying  change  in  the  flowfield.  However,  the  2-D 
computations  and  experiments  have  demonstrated  both  solutions[30,  18]. 

Figure  1.9  is  the  schematic  of  Schmisseur  and  Gaitonde’s  2-D  solution  domain[30]. 
The  2-D  slice  in  Fig.  1.4  is  cut  in  half  with  a  symmetry  boundary  condition  down  the 
center  to  reduce  the  number  of  cells  and  corresponding  computation  time.  Their  com¬ 
putational  solutions  employed  compression  angles  of  22°  <  9  <  28°  at  a  freestream 
Mqo  =  4.96  and  Reynolds  number  based  on  w  of  2  million.  These  conditions  corre¬ 
spond  to  those  at  which  experimental  data  were  obtained  by  Ivanov,  et  al.  [17,  20] 
in  order  to  facilitate  validation  of  both  Navier-Stokes  and  Euler  simulations. 


Figure  1.9  Schematic  of  Schmisseur  and  Gaitonde’s  2-D  physical  domain[30] 

The  hysteresis  phenomenon  shown  in  Fig.  1.10  was  simulated  by  Schmisseur 

and  Gaitonde  for  inlet  conditions  with  g/w  =  0.42  and  g/w  =  0.34[30].  The  half 

throat  diameter  is  defined  g  and  w  is  the  length  of  the  compression  surface  as  shown 

in  Fig.  1.9.  For  the  g/w  =  0.34  results  in  the  dual-solution  domain,  the  compression 

angle  was  incrementally  increased  from  the  theoretical  von  Neumann  condition,  6^  = 

20.87°.  The  transition  from  RR  to  MR  occurred  at  a  compression  angle  between 
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27.5°  and  27.85°.  This  angle  range  is  in  close  agreement  with  the  theoretical  value 
for  detachment,  9 d  =  27.7°.  Once  a  MR  solution  was  obtained,  the  solution  was 
used  as  the  initial  condition  for  solution  of  the  lower  compression  angle.  The  MR 
pattern  persisted  through  much  of  the  dual-solution  space  transitioning  to  an  RR 
configuration  before  the  von  Neumann  condition  was  reached.  The  transition  from 
a  MR  to  RR  solution  occurred  between  9  =  22°  and  23°.  The  inviscid  numerical 
results  for  the  height  of  the  normal  shock  Mach  stem  of  Schmisseur  and  Gaitonde[30] 
showed  good  agreement  with  the  experimental  data  and  Euler  solutions  by  Ivanov,  et 
al.  [20]  as  shown  in  Fig.  1.11.  This  figure  indicates  particularly  good  agreement  in  the 
region  of  the  dual-solution  domain  near  detachment.  In  this  region  the  experimental 
data  generally  lie  within  the  mesh  resolution  of  the  numerical  values. 


Increasing  Compression  Ramp  Angle 


Decreasing  Compression  Ramp  Angle 
Increasing  Compression  Ramp  Angle 


Decreasing  Compression  Ramp  Angle 


Figure  1.10  Hysteresis  loop  for  Mach  5,  2-D  solutions[30] 
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Shock  Wave  Angle  (Degrees) 

Figure  1.11  Variation  of  Mach-stem  height  with  6  for  2-D  solutions [30]  (Medium 
corresponds  to  a  grid  resolution  of  206  —  329  x  137) 

Ivanov,  et  al.  produced  the  hysteresis  experimentally  for  a  freestream  Mach  of 
4 [?] .  The  results  demonstrated  the  dual-solution  domain  from  6  =  34°  to  0  —  38.2°. 
These  values  are  very  close  to  the  von  Neumann  and  detachment  angles  for  =  4 
of  33.4°  and  39.2°,  respectively. 

Since  the  transition  from  regular  reflection  to  Mach  reflection  can  initiate  many 
problems,  control  and/or  prevention  of  this  flowfield  transition  is  an  obvious  goal. 
Researchers  have  examined  the  factors  which  transition  the  flow  and  methods  of 
controlling  the  transition [13,  14,  19,  34,  35].  The  experiments  of  Yan,  et  al.  with 
laser  energy  deposition  have  successfully  transitioned  a  RR  to  MR  and  a  MR  to 
RR[34,  35].  This  method  energizes  the  flow  in  front  of  the  shock  system.  In  the 
current  research,  magnetogasdynamics  is  investigated  as  a  method  of  controlling 
this  flow  transition. 
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1.3  Magnetogas  dynamics 

Magnetofluidmechanics  (MFD)  unifies  electrodynamics,  thermodynamics,  and 
fluid  mechanics  into  one  discipline,  the  study  of  the  interaction  between  electro¬ 
magnetic  fields  and  a  continuum  fluid [15].  These  interactions  can  be  modelled  by 
coupling  the  Maxwell  equations  with  the  Navier-Stokes  equations.  Magnetogasdy- 
namics  (MGD)  is  the  subset  of  MFD  associated  with  continuum  gases  and  the  many 
environmental  variables  present  in  hypersonic  flight.  Applying  MFD  to  the  hyper¬ 
sonic  flight  regime  allows  many  simplifying  assumptions  to  be  made.  The  induced 
magnetic  field  is  negligible  for  low  magnetic  Reynolds  number  flows.  The  flow  is 
weakly  ionized  so  the  conductivity  is  low.  The  low  conductivity  specifies  a  mag¬ 
netic  Reynolds  number  much  less  than  1.  This  simplification  allows  the  use  of  the 
low  magnetic  Reynolds  number  form  of  the  governing  equations.  The  MGD  source 
terms  that  remain  are  the  ponderomotive  forces  and  accompanying  energy  interac¬ 
tion  terms  [15]. 

In  recent  years,  magnetogasdynamics  has  been  suggested  as  a  flow  control 
mechanism  for  the  hypersonic  flow  regime[5,  7,  8,  22,  23,  25].  Researchers  have 
examined  the  use  of  electromagnetic  fields  to  control  shock  location  for  inlet  mass- 
capture,  enhance  fuel-air  mixing,  suppress  boundary  layer  transition,  and  serve  as  a 
means  to  efficiently  withdraw  energy  from  the  flow  in  the  inlet  and  add  energy  to 
the  flow  in  the  nozzle  for  thrust  enhancement.  The  primary  mechanisms  are  pon¬ 
deromotive,  or  Lorentz,  forces  and  accompanying  energy  interactions,  which  include 
both  reversible  and  irreversible  components.  The  potential  use  of  MGD  to  transition 
a  SCRAMjet  inlet  from  Mach  reflection  to  regular  reflection  is  an  enticing  prospect. 
In  theory,  a  Lorentz  force  applied  tangent  to  the  flow  could  be  used  to  reduce  or  even 
eliminate  the  MR.  This  application  of  the  Lorentz  force  reduces  the  normal  shock 
stem  height,  helps  to  sustain  supersonic  combustion,  and  minimizes  the  chance  of 
engine  unstart.  The  same  applied  force  could  also  prevent  transition  to  MR  before 
it  could  occur. 
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1- 4  Current  Research 

The  current  research  effort  begins  with  the  regular  and  Mach  reflection  data 
of  Schmisseur  and  Gaitonde[30]  and  applies  magnetogasdynamic  forces  to  determine 
the  sensitivity  of  both  types  of  reflections  to  decelerating  and  accelerating  Lorentz 
forces.  A  2-D  inviscid  inlet  is  examined  at  multiple  compression  ramp  angles  in  the 
dual-solution  domain.  A  converged  RR  is  subjected  to  a  decelerating  Lorentz  force 
to  determine  its  sensitivity  to  transition  to  a  MR  and  possible  unstart.  In  addition, 
an  accelerating  Lorentz  force  is  applied  to  a  MR  to  try  and  transition  the  MR  to  a 
more  efficient  RR. 

This  numerical  study  employs  a  3-D  magnetogasdynamic  Computational  Fluid 
Dynamics  (CFD)  code [3,  4,  5,  6,  7,  8].  The  code  was  written  by  the  Computational 
Branch  of  the  Aerodynamics  Division  of  the  Air  Vehicles  Directorate  (VAAC)  within 
the  Air  Force  Research  Laboratory  (AFRL).  The  code  solves  the  Euler  or  Navier- 
Stokes  equations  with  the  addition  of  source  terms  for  the  MGD  forces  and  energy 
interactions.  The  AFRL /VAAC  code  has  been  verified  for  various  problems  including 

2- D  flow  past  a  flat  plate  and  a  laminar  MGD  boundary  layer  flow[5,  7].  For  the 
current  research,  the  inviscid  Euler  equations  are  solved. 

The  magnetogasdynamic  governing  equations  used  in  this  research  are  devel¬ 
oped  in  Section  2  from  the  Navier-Stokes  equations  and  the  Maxwell  equations.  In 
section  3,  the  numerical  algorithms  used  to  solve  the  magnetogasdynamic  equations 
are  presented.  Section  3  also  includes  a  detailed  explanation  of  the  the  computational 
domain  set-up  for  this  research.  Section  4  begins  with  duplication  of  Schmisseur  and 
Gaitonde’s  hysteresis.  The  sensitivity  of  the  regular  reflection  to  a  decelerating 
Lorentz  force  and  the  Mach  reflection  to  an  accelerating  Lorentz  force  are  examined, 
including  the  dependency  on  location  of  the  applied  field  and  the  magnitude  of  the 
applied  magnetic  and  electric  fields.  Section  5  presents  some  conclusions  on  the  use 
of  MGD  to  control  Mach  reflections  and  possible  future  work. 
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2.  Governing  Equations 


2.1  Navier-Stokes  Equations 

The  Navier-Stokes  equations  (NSE)  are  the  governing  equations  of  fluid  flow. 
These  equations  are  based  on  the  laws  of  the  conservation  of  mass,  momentum,  and 
energy.  When  radiative  heat  transfer  and  body  forces  are  neglected,  these  laws  can 
be  expressed  as  the  continuity,  momentum,  and  energy  equations  below: 


d_P_  ,  d(pbj) 

dt  dxi 


(2.1) 


(2.7) 


Kf  «  Kq,/ 


T  \  5  T0  +  5, 


Tn 


M/ 


T  +  5K/ 


where  /i0,/  =  1.716  x  10  ,  T0  =  273  K,  S)(/  =  111  K,  /c0,/  =  2.41  x  10  ,  and 
SKf  =  194  K  for  air [33].  The  total  enthalpy  and  total  energy  which  appear  in  the 
NSE  are  defined  as  follows: 


ht  =  h+^VM  (2.8) 

e*  =  ht~-  (2.9) 

2.1.1  Ideal  Gas  Law.  The  NSE  as  represented  above  have  six  unknowns, 
(p,u,v,w,p,et),  and  5  equations.  Assuming  that  air  is  a  calorically  perfect  gas  adds 
the  sixth  equation,  the  ideal  gas  law. 


p  =  plZT  (2.10) 

In  the  preceding  equation,  7 Z  is  the  universal  gas  constant  for  air. 

2.1.2  Vector  Form  of  the  Navier-Stokes  Equations.  The  NSE  can  be  writ¬ 
ten  in  conservative  vector  form  as  follows: 


dU  dE  dE  dG 
dt  dx  dy  dz 

where  U  is  the  vector  of  conservative  variables: 


(2.11) 


U  —  [p  pu  pv  pw  pet\ 


T 


(2.12) 
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E.  F.  and  G  are  the  total  fluxes  in  the  x,  y,  and  z-directions,  respectively,  which  can 
be  broken  down  into  inviscid  and  viscous  flux  vectors: 


E  —  Ei  —  Ev 


F  =  Fi-Fv 
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(2.13) 


(2.14) 


(2.15) 


and  S  is  the  source  term: 


S  = 


0 

0 

0 

0 

0 


(2.16) 


The  inviscid  Euler  equations  are  formed  by  neglecting  the  viscous  fluxes. 
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2.2  Maxwell  Equations 

The  governing  equations  of  electromagnetics  are  known  as  the  Maxwell  equa¬ 
tions  (ME).  These  equations  are  valid  in  all  reference  frames,  whether  observation 
is  in  the  laboratory  frame  with  the  media  moving  relative  to  the  observer,  or  the 
rest  frame  where  the  observer  is  at  rest  compared  to  the  media,  denoted  by  a  prime. 
These  equations  are  defined  as  follows  and  are  changed  to  the  rest  frame  by  the 
addition  of  a  prime  to  all  terms[15]: 


v  •  E  =  — 

(2.17) 

ee 

V-B  =  0 

(2.18) 

d  B 

(2.19) 

VxE  =  ~m 

d  E 

(2.20) 

VxB  =  Hej  +  /Teee  — 

where  E  is  the  electric  field,  B  is  the  magnetic  flux  density,  j  is  the  conduction 
current  density,  /xe  is  the  permeability,  and  ee  is  the  permittivity. 

2.2.1  Ohm’s  Law.  Ohm’s  law  is  defined  for  linear  isotropic  media  in  the 
rest  frame  as  follows: 


j'i  =  Vji&j 


(2.21) 


where  a  is  the  conductivity. 

The  Maxwell- Lorentz  transformations  (MLT)  are  used  to  transform  the  Maxwell 
equations  from  the  rest  frame  to  the  laboratory  frame.  Several  terms  are  added  to 
Ohm’s  law  during  the  transformation  to  allow  for  a  greater  range  of  media.  When 
the  media  is  nonuniform  in  all  directions,  anisotropic,  and  V2  <<  c2,  Ohm’s  law  for 
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a  nonlinear  media  is  as  follows  [15]: 


ji  —  &ji[Ej  +  (V  x  B)j]  +  peVi  (2.22) 

where  pe  is  the  space  charge. 

2.2.2  Constitutive  Relations.  The  constitutive  relations  for  linear  isotropic 
dielectrics  and  magnetic  materials  expressed  in  the  rest  frame  are  as  follows[15]: 


(2.23) 

ee 

B’  =  peH' 

(2.24) 

where  D  is  the  electric  flux  density  and  H  is  the  magnetic  field. 

The  MLT  are  used  to  transform  the  constitutive  relations  to  the  laboratory 
frame.  When  V2  <<  c2  and  the  medium  is  isotropic,  the  constitutive  relations  are 
as  follows  [15]: 


D  —  ee[E  +  (1 
B  =  pe[H-(  1 


1 


€e 

Me 

€0,e 

VO,e 

1 

_€e_ 

Me 

£0,e  MO,e 


)V  x  B] 
)V  x  D\ 


(2.25) 

(2.26) 


2.2.3  Magnetogas  dynamic  Assumptions.  When  an  electrically  conducting 
fluid  is  moving  in  the  presence  of  a  magnetic  field,  the  flow  of  the  fluid  is  influenced 
by  the  field  and  the  field  is  influenced  by  the  moving  fluid.  If  certain  assumptions 
are  made,  these  interactions  can  be  expressed  as  equations  which  govern  magnetoflu- 
idmechanic  flow.  When  these  assumptions  are  applied  to  a  gas  such  as  air,  we  apply 
the  scientific  name  magnetogasdynamics  and  these  assumptions  become  the  MGD 
assumptions.  The  MGD  assumptions  are  as  follows  [15]: 
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MGD  Assumption  1:  |V"|2  <<  c2,  the  magnitude  of  the  velocities  dealt  with  in 
fluid  dynamics  are  much  less  than  the  speed  of  light  allowing  the  1  —  ( V/c 2) 
term  to  be  set  to  unity. 

MGD  Assumption  2:  E  ~  0(V  x  B ),  the  electric  held  is  of  the  order  of  any 
induced  effects  which  implies  that  the  applied  magnetic  held  is  much  greater 
than  the  induced  magnetic  held. 

MGD  Assumption  3:  ^  ~  0,  ^  is  also  zero  by  the  constitutive  relations.  This 
assumption  then  implies  that  V  x  B  =  /iej. 

MGD  Assumption  4:  eE 2  <<  — ,  the  electric  energy  is  insignificant  compared  to 
the  magnetic  energy. 

MGD  Assumption  5:  j  =  cr(E  +  V  x  B ),  the  conductivity  is  also  considered 
independent  of  magnetic  held  and  constant  with  frequency.  This  assumption 
implies  that  j'  =  j. 

MGD  Assumption  6:  /  =  peE  +  j  x  B. 

2.2.4  Maxwell  Equations  for  Magnetogasdynamic  Flow.  When  the  ap¬ 
proximations  above  are  applied  to  the  ME,  a  new  set  of  equations  which  govern 
electromagnetic  helds  are  formed.  These  equations  in  the  laboratory  frame  are  as 
follows  [15]: 


(2.27) 

V  X  B  =  p.ej 

(2.28) 

V-J  =0 

(2.29) 

V-B  =  0 

(2.30) 

j  =  a(E  +  V  x  B) 
2-6 


And,  Ohm’s  law  becomes: 


(2.31) 


All  subsequent  equations  are  written  in  the  laboratory  frame  unless  otherwise 


noted. 


2.3  Magnetogas  dynamic  Equations 

The  Magnetogasdynamic  Equations  (MGDE)  result  from  the  combination  of 
the  Maxwell  equations  for  MGD  and  the  NSE.  These  equations  describe  the  interac¬ 
tion  between  electromagnetic  fields  and  electrically  conducting  gases  in  a  continuum 
governed  by  the  MGD  assumptions.  The  MGDE  are  defined  as  follows: 


dp  d(pVj) 
dt  dxi 


0 


d(pVj)  d^pVjVj)  _  d(p5ij )  _  dr^ 

dt  dxj  dxj  dxj 

d(pet)  d(phtVj)  _  d(gj  +  TjjVf) 

dt  dxj  dxj 


jxB 


E-j 


(2.32) 

(2.33) 

(2.34) 


where  (j  x  B)  is  the  Lorentz  force,  and  E-j  is  the  resulting  energy  interaction. 


2.3.1  Non-dimensionalization  of  the  Magnetogasdynamic  Equations.  By 
non-dimensionalizing  the  equations,  problems  of  different  magnitude  scales  can  be 
compared  independent  of  the  reference  conditions.  The  dimensionless  quantities 
used  in  the  normalization  of  the  MGDE  are  as  follows  where  an  asterisk  represents 
a  dimensionless  quantity [15]: 


L* 

P* 

E* 


L 

V 

Lo 

P 

V* 

V0 

T 

Po 

rj~\ * 

E 

B 

Eo 

B* 

B~o 

(2.35) 
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These  initial  dimensional  quantities  of  L0,  Vo,  po,  T0,  E0,  and  B0  are  used 
to  non-dimensionalize  the  other  variables,  and  all  dimensionless  variables  are  sub¬ 
stituted  into  the  MGDE.  Several  dimensionless  parameters  are  formed  during  the 
non-dimensionalization  which  help  to  characterize  the  MGD  problem.  The  Reynolds 
number  (Re)  is  the  ratio  of  inertial  forces  to  viscous  forces,  defined  as  follows: 

Re  =  M £  (2.36) 

uf 

The  magnetic  Reynolds  number  (Rm)  is  a  measure  of  the  magnitude  of  the  induced 
magnetic  field  compared  to  the  total  magnetic  field: 


Rm  =  V0L0a0po,e  (2.37) 

The  Mach  number  (M)  and  magnetic  Mach  number  (Mm)  are  defined  as  follows: 


o  o 


Mm  = 


boy /Po7hp= 


(2.38) 

(2.39) 


The  Prandtl  number  (Pr)  is  the  ratio  of  kinematic  viscosity  to  thermal  diffusivity: 


Pr  =  (2.40) 

K, 

The  magnetic  Prandtl  number  (Pm)  is  the  ratio  of  vorticity  diffusion  to  magnetic 
diffusion: 


Pm  =  (7{)vfiM),e 


Rm 

Re 


(2.41) 
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The  interaction  parameter,  Q,  is  used  when  Rm  «  1  and  is  defined  as  the  ratio 
of  the  ponderomotive  force  to  the  inertial  force  and  is  0(1).  By  increasing  Q  the 
magnitude  of  the  magnetic  field  is  increased. 


<j0BqL0  Rm 
PoVo  M?n 


(2.42) 


The  load  factor,  /C,  is  defined  as  the  ratio  of  the  electric  field  to  the  induced  effects. 
An  increase  in  /C  signifies  a  larger  applied  electric  field,  Eq. 


1C 


-Eq 

B0V0 


(2.43) 


The  non-dimensional  MGDE  are  identical  to  their  dimensional  counterparts 
except  for  the  use  of  non-dimensional  quantities  and  the  differences  in  the  stress 
tensor,  heat  flux,  and  MGD  source  terms. 


* 

Re 


-  dvi  dvi  2 5  d_vi 

_  dx*  dx*  3  v  dx*k  _ 

p)  dT* 

(7  —  1  )M%0RePr  dx* 


(2.44) 

(2.45) 


The  MGD  source  terms  are  now  multiplied  by  the  interaction  parameter.  The  re¬ 
sulting  non-dimensional  Lorentz  force  is  Q(j*  x  B*)  and  the  energy  interaction  is 
QE*  ■  j* .  The  non-dimensional  perfect  gas  law  is  defined: 


with 


7  M^oP* 

p* 


(2.46) 


P*  =  (7  -  1  )p*e* 


(2.47) 
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For  the  remainder  of  this  document  the  *  is  dropped  and  all  quantities  are  dimen¬ 
sionless  unless  otherwise  noted. 

2.3.2  Vector  Form  of  the  Magnetogasdynamic  Equations.  The  equations 
resulting  from  the  non-dimensionalization  can  then  be  written  in  a  conservative 
vector  form  with  the  electromagnetic  effects  in  the  source  term  as  follows: 


dU  DE  <9F  dG  _ 
dt  dx  dy  dz 


(2.48) 


The  only  change  from  the  non-dimensional  conservative  vector  form  of  the  NSE  is 
in  the  source  term  S,  where: 


0 

Q(jyBz  -  jzBy) 
Q(jzBx  -  jxBz) 

Q^jxBy  jy  Bx  ) 

Q(EXjX  T  Fyjy  T  Ezjz) 


(2.49) 


2.3.3  Curvilinear  Transformation.  In  order  to  accurately  model  most  phys¬ 
ical  problems,  the  grid  created  for  the  physical  domain  has  non-uniform  spacing  in  at 
least  one  direction.  This  mesh  might  also  be  body-fitted  involving  a  non-orthogonal 
coordinate  system.  To  account  for  this,  the  vector  form  of  the  MGDE  presented 
above  is  transformed  to  a  generalized  coordinate  system  (£,t/,C)- 


f  =  £{x,y,z) 

V  =  r}{x,y,z)  (2.50) 

C  =  C  (x,y,z) 

T  =  t 
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The  differential  expressions  for  this  coordinate  transform  are  as  follows: 


=  ixdx  +  tydy  +  izdz 

dr]  =  r]xdx  +  r]yd/y  +  rjzdz  (2-51) 

d(  =  C Xdx  +  Cydy  +  Czdz 

dr  =  dt 

Similarly,  the  differential  expressions  for  the  inverse  transformation  are  as  follows: 


dx  =  x^d^  +  xvdr]  +  x^d(  +  xTdr 

dy  =  yxidi  +  yvdr /  +  ycd(  +  yTdr  (2.52) 

dz  =  zxid£  +  Zffdr /  +  z^d(  +  zTdr 
dt  =  dx 

The  transformation  metrics  are  defined  using  Equation  2.51  and  Equation  2.52[31]: 


=  J  (VvzC  ~Vczv) 

C,y  ./  (XrjZ^  X^Z'q')  (2.53) 

=  J (l/riVC  ~  yCVri) 

Vx  =  -J(y£ZC-y{Zz) 

r]y  =  J(x%Z£  —  x^zg)  (2-54) 

Vx  =  -J(ym~y<:yd 
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Ca;  Ur 


C y  =  ~  J(x^Zr,  ~  XvZz)  (2.55) 

Cx  J  (y^yv  yvyd 


where  the  Jacobian  of  transformation,  J ,  is  defined  as  follows [31]: 


J  = 


x^(yvz c  -  vczv)  -  xv(ytz c  -  vch)  +  xdyzzv  -  y^n) 


(2.56) 


Applying  the  transformations  to  Equation  2.48  and  manipulating  the  equation 
into  strong  conservation  form,  the  new  terms  U,  E,  F,  G,  and  S  are  formed. 


E  =  I(E^  +  F4  +  G^) 

F  =  —  {Eir]x  +  Frjy  +  Grjz)  (2-57) 

G  =  j(E(x  +  FCy  +  GCz) 


with  the  resulting  vector  equation[31]: 


dU  <9E  5F  dG 
dr  <9£  drj  d( 


(2.58) 
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3.  Numerical  Implementation 

This  section  presents  the  discretization  of  the  Magnetogasdynamic  Equations  as  they 
are  employed  in  the  AFRL/VAAC  three-dimensional  MGD  CFD  code[3,  4,  5,  6,  7,  8]. 
For  this  inviscid  research  study,  the  viscous  terms  in  the  MGDE  are  neglected.  The 
curvilinear  vector  form  of  the  MGDE,  Eqn.  2.58,  on  page  2-12  is  rewritten  as  follows: 


dU  dE  dF  dG 
dr  dr/  d( 


where  the  fluxes  are  inviscid 


(3.1) 


E  =  jiE^  +  FNy  +  CnQ 

F  =  j(EIrjx  + F^y  +  G^)  (3.2) 

G  =  —  (Ej(x  +  Fj(y  +  GNz) 

Equation  3.1  is  solved  using  a  nodal  finite- volume  approach.  The  physical  grid 
is  transformed  to  the  evenly  spaced  computational  grid  using  Eqns.  2.53,  2.54,  and 
2.55.  The  discretization  of  Equation  3.1  has  the  form: 

dU 

—  =  -K(U)h,k  (3-3) 

where  R  is  the  residual  vector: 


T>  (TT\n  _  T?71  _i_  fpn 

)i,j,k  1/2, j,k  i,j+l/2,k  1/2, k 


+G' 


i,j,k-\- 1/2 


G 


i,j,k— 1/2 


+« 


(3.4) 


ij,k 
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The  indices  (i,j,  k )  represent  the  node  values,  and  the  indices  (i  ±  1/2,  j,  k ),  ( i,j  ± 
1/2,  k),  and  (i,j,k  ±  1/2)  represent  the  right/left,  top/bottom,  and  front/back  cell 
interfaces,  respectively. 

3.1  Flux  Representation 

The  fluxes  across  the  boundaries  of  each  cell,  E,  F,  and  G,  are  defined  as 
follows[31]: 


E  = 


1 

J 


pU 

puU  +  £xp 

PVU  +  £yP 
pwU  +  £zp 
phiU 


F  = 


1 

J 


PV 

puV  +  T]xp 
pvV  +  rjyP 
pwV  +  r]zp 
phtV 


G 


1 

J 


pW 

puW  +  CxP 
pi’W  +  CyP 
pwW  +  (zp 
phtW 


(3.5) 


(3.6) 


(3.7) 
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where  U,  V,  and  W  are  the  contravariant  velocity  components  in  the  £,  77,  and 
directions,  respectively,  dehned  as  follows: 


U 

V 

W 


fxF  +  iyV  F  f,zW 

(3.8) 

VxU  F  rjyV  F  rjzw 

C xU  F  C yV  F  (ZW 

(3.9) 

3.1.1  Roe  Flux- Difference  Splitting  Scheme.  The  flux  vectors  are  dis¬ 
cretized  using  Roe’s  upwind  method  with  a  van  Leer  harmonic  limited  Monotone 
Upstream-centered  Schemes  for  Conservation  Laws  (MLISCL)  variable  extrapolation [28, 
32,  11],  The  Roe  scheme  is  a  flux-difference  splitting  method  created  to  more  ac¬ 
curately  capture  the  discontinuities  present  in  high-speed  fluids.  Roe’s  algorithm  is 
applied  to  each  coordinate  independently.  The  following  derivation  applies  the  Roe 
scheme  in  the  £  direction.  The  cell  interface  flux  is  written  as  the  exact  solution  to 
an  approximate  Riemann  problem  [32],  given  by 


Ei+l/2,j,k 


2  L 


E{Ul)  +  E{Ur)-\A\{Ur-Ul) 


i+l/2,j,k 


where  A  is  the  flux  jacobian 


(3.10) 


A  = 


d  E 
dU 


The  Roe  vector  in  Eqn.  3.10  can  now  be  defined  as  follows [32]: 


(3.11) 


A|(C/„  -  UL)  =  |A|AC7  = 


OL  4 

ua4  F  kxa§  T  o:g 
vex 4  T  kyCx^  T  07 
Wa4  T  kza§  T  cxg 

~  ~  ~  _  a2 

Ha4F  uoi§  F  ua§  F  vcxi  F  was - ai 
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where  the  a  terms  are  defined 
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(3.13) 


The  interface  areas  and  normal  velocity  terms  are  defined  as  follows: 

7  Sir 

**  “  Ml 

P  - 

y  ~  |vei 

7,  _ 

'  |V{| 

u  =  kxu  +  +  kzw 


(3.14) 
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In  the  above  equations,  the  Roe  averaged  variables  are  defined  as  follows: 


Ri+ 1/2 
Pi+ 1/2 
fij+1/2 

fij+1/2 

Vi+l/2 

Wj+1/2 

Hi+ 1/2 
~2 

ai+l/2 


[PR 

V  PL 

PlRi+ 1/2 

+  unRi+i/2 

1  +  Ri+\/2 
Pl  +  fiR-Rj+1/2 
1  +  Ri+l/2 

VI  +  hR-Ri+1/2 

1  +  Ri+l/2 

WL  +  Wfi-Rj+1/2 

1  +  -Rj+1/2 
Hl  +  ^R-Rj+1/2 
1  +  -Rj+l/2 

(7  -  l)[#j+i/2  -  (fi2  +  fi2  +  w2)/2] 


(3.15) 


The  Roe  scheme  is  capable  of  producing  a  non-physical  expansion  shock  given 
certain  conditions.  The  expansion  shock  violates  the  second  law  of  thermodynamics. 
In  order  to  prevent  this  occurrence,  an  entropy  fix  is  added  to  the  Roe  scheme  [31]. 
For  the  problem  examined  in  the  current  research,  the  conditions  needed  to  form  an 
expansion  shock  are  not  present. 

3.1.2  MUSCL  Higher-Order  Variable  Extrapolation.  In  order  to  achieve  a 
higher  order  of  spatial  accuracy,  the  primitive  variable  vector  to  the  left  and  right  of 
the  cell  interface  is  extrapolated.  The  primitive  variable  vector,  V,  is  defined 


V 


P 

u 

v 

w 


V 


(3.16) 
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The  higher  order  left  and  right  cell  interface  values  of  the  primitive  variables 
are  calculated  using  a  van  Leer  harmonic  limited  MLTSCL  variable  extrapolation [11], 
The  limiter  prevents  the  slope  across  the  interface  from  exceeding  stability  bounds. 


where 
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(3.17) 


(3.18) 


(3.19) 


(3.20) 
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In  Eqns.  3.17  and  3.18,  /q  is  set  to  1/3  for  the  higher-order  calculations  and  the 
left  and  right  cell  interface  values  are  set  to  the  i  and  i  +  1  node  values,  respectively, 
for  first-order  calculations. 

3.2  Temporal  Discretization 

The  code  includes  an  explicit  Runge-Kutta  time  integration  method  and  an 
implicit  Beam- Warming  time  integration  method.  Explicit  schemes  severely  limit 
the  time  step-size  based  on  the  size  of  the  cells  in  the  physical  domain.  These 
limitations  are  made  worse  near  the  wall  where  the  mesh  is  very  hne  and  there 
are  large  stretch  factors[6].  To  avoid  this,  an  implicit  method  with  a  sub-iteration 
strategy  is  incorporated  to  solve  the  low  magnetic  Reynolds  number  form  of  the 
governing  equations.  This  integration  is  covered  in  detail  by  Gaitonde  and  Poggie  in 
Ref.  [6],  including  validation  of  the  implicit  scheme.  The  basic  formula  to  compute 
the  change  in  the  solution  vector,  A U,  is: 
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1  +  0 
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J-x  + 
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j-x  + 
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A  U  = 


J-i(l+^)Vi’-(l  +  2^)V"  +  ^A'-1  + 
At 


(3.21) 


In  Eqn.  3.21  At  is  the  generalized  time-step  size,  5%,  5V,  and  5$  denote  central 
second-order  difference  operators  and  0  controls  the  temporal  order  of  accuracy  of 
the  scheme.  A,  B  and  C  are  transformed  flux  Jacobians  and  R  is  the  residual 
evaluated  with  the  Roe  scheme. 

Sub-iterations,  indicated  by  the  counter  p  in  Eqn.  3.21,  reduce  linearization, 
factorization,  and  explicit  boundary  condition  implementation  errors.  Each  iteration 
typically  requires  two  or  three  sub- iterations  to  achieve  inner-iteration  convergence  [6] . 
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3.3  Computational  Set-up 

Figure  3.1  illustrates  the  physical  domain  used  for  this  work.  Increases  in  6 
were  accomplished  by  moving  the  leading  edge  of  the  ramp  down  and  back,  effectively 
opening  up  the  inlet.  The  ratio  g/w  was  held  constant  at  0.34  by  allowing  the  ramp 
to  be  hinged  at  the  throat.  Fine  213  x  99  grids  were  created  for  all  compression 
angles,  22°  <  9  <  28°. 


Figure  3.1  Schematic  of  the  two-dimensional  domain 


The  duct  angle  in  Fig.  3.1  is  0°.  This  difference  from  Schmisseur  and  Gaitonde’s 
domain  [30]  shown  in  Fig.  1.9  is  due  to  inaccuracies  encountered  with  a  duct  angle  of 
25.4°.  The  expansion  angle  created  by  combining  6  and  the  duct  angle  of  25.4°  was 
too  large  for  inviscid  flow.  This  expansion  angle  caused  the  computational  code  using 
Roe  flux-differencing  to  become  unstable.  The  early  work  with  the  25.4°  duct  angle 
used  the  finite-compact-differencing  scheme  with  dampening  presented  by  Gaitonde 
in  Ref.  [3].  This  numerical  scheme  introduced  artificial  viscosity  for  dampening, 
allowing  the  inviscid  flowfield  to  converge  inaccurately.  Discovery  of  this  problem 
led  to  the  0°  duct  angle  shown  in  Fig.  3.1. 

The  computations  for  this  inviscid  study  were  performed  at  a  Mach  number 

of  4.96.  In  order  to  completely  determine  the  freestream  variables,  the  Reynolds 

number  based  on  w  was  set  to  2  million.  A  coarse,  fine,  and  ultra-fine  mesh  were 
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created  for  the  6  =  26°  ramp  angle.  These  grids  were  used  for  the  mesh  refinement 
study.  Figure  3.2  shows  the  coarse,  fine,  and  ultra-fine  26°  mesh.  The  line  running 
near  the  top  of  the  domain  is  the  symmetry  plane.  The  coarse  grid  is  107  x  50,  and 
the  ultra-fine  grid  is  425  x  197.  Each  grid  refinement  step-up  corresponds  to  exactly 
twice  the  cells  in  both  the  £  and  rj  directions.  All  grids  include  a  mirror  image  of 
the  top  four  rows  about  y  =  0.  These  extra  rows  are  used  for  a  symmetry  boundary 
condition.  The  computations  included  5  total  planes  in  £  direction.  The  two  planes 
on  each  side  of  the  center  plane  are  copies  of  the  center  plane  forcing  the  solutions 
to  be  two-dimensional. 

The  magnetogasdynamic  equations  resolve  the  flow  within  the  domain,  but 
boundary  conditions  must  be  specified  to  define  the  flow  conditions  along  the  com¬ 
putational  boundary.  The  MGD  code  reverts  to  first  order  at  the  boundaries.  This 
facilitates  solution  of  a  domain  without  resorting  to  non-physical  ’’ghost”  nodes  to 
create  higher-order  boundary  conditions. 

For  the  current  research,  the  entrance  of  the  inlet  is  an  inflow  condition  and 
the  exit  is  an  outflow  condition.  The  walls  of  the  inlet  are  specified  as  an  inviscid-slip 
boundary  condition  and  the  inlet  is  symmetric,  so  only  half  the  inlet  is  modelled  to 
decrease  the  number  of  nodes  and  corresponding  computational  time.  A  symmetry 
boundary  condition  (BC)  is  applied  along  the  split. 

The  hypersonic  flow  conditions  allow  the  use  of  a  supersonic  inflow  BC  on 
the  entrance  to  the  inlet.  The  values  of  the  primitive  variables  on  the  boundary, 
V(i,  are  equal  to  the  freestream  values  of  p,  u,  v,  w,  and  p.  With  the  non- 
dimensionalization,  p  =  u  =  1,  v  =  w  =  0,  and  p  =  Po/(poWq)-  Figure  3.3  illustrates 
the  supersonic  inflow  BC. 

The  flow  exiting  the  inlet  should  remain  supersonic  due  to  the  expansion  after 
the  compression  ramp.  This  allows  specification  of  a  supersonic  outflow  BC  at  the 
exit  of  the  inlet.  The  values  of  the  primitive  variables  on  the  boundary,  V (ni,j.k)i  are 

equal  to  the  primitive  variable  values  one  node  in,  ~V(ni-i,j,k),  as  shown  in  Fig.  3.3. 
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Figure  3.2  26°  coarse,  fine,  and  ultra-fine  grids  with  corresponding  symmetry  line, 

g/w  =  0.34 
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Outflow 

Boundary 


Figure  3.3  Schematic  of  the  supersonic  inflow  and  outflow  BCs 


The  values  are  simply  copied  from  the  (m  —  1,  j,  k )  location  to  the  boundary  location 
at  ( ni,j,k ). 

The  symmetry  boundary  condition  simply  specifies  that  the  nodes  above  the 
symmetry  plane  are  mirrored  from  the  nodes  below  the  plane.  If  the  symmetry  plane 
falls  at  ( i,nj,k ),  the  primitive  variables  one  node  above,  V^nj+ i,k),  are  set  to  the 
values  one  node  below,  V  (i,nj-  1>k).  The  trend,  V ^nj+m,k)  =  V (i,nj-m,k),  continues 
incrementing  until  m  =  4.  The  symmetry  BC  at  the  center  of  the  inlet  decreases  the 
number  of  nodes  by  almost  50%,  correspondingly  decreasing  the  computation  time 
for  convergence. 

The  inviscid  slip  BC  along  the  surface  of  the  inlet  projects  the  velocity  vector 
one  node  from  the  wall,  U^:2,k)  onto  the  surface  of  the  wall,  as  shown  in  Fig.  3.4. 
The  velocity  in  the  z-direction,  w,  is  set  to  zero.  This  process  starts  with  calculation 
of  the  angle  of  the  wall,  ©#,  as  follows 


Qb  =  arctan  (  V{i  1,1,fc) — }!&h!*L\  (3.22) 

\  *£(i— 1,1, fc)  £(z,l,A0  / 

The  velocity  vector,  U^k)  at  one  node  in  is  then  projected  onto  the  wall.  The 
magnitude  of  this  velocity  is  calculated  as  follows 


||^(vL,fc)||  =  cos(0B)-u(ii2,fc)  +  sin(0B)up>2,fc) 
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(3.23) 


The  u  and  v  components  of  U^t i  k)  shown  in  Fig.  3.4  are  then  calculated 


U(i,i,k)  =  cos(0B)||f/(i.1,fc)||  (3.24) 

V(i,i.k)  =  sin(0B)||£/(j)i!fe)|| 


Figure  3.4  Schematic  of  the  inviscid  slip  BCs 


The  pressure  and  density  at  the  specified  wall  location  are  set  equal  to  the 
pressure  and  density  one  node  in,  P(i,i,k)  =  P(i,2,k)  and  P(i,i,k)  =  P(i,2,k)- 

For  all  compression  angles  9  <  28°,  the  regular  symmetry  solutions  were  ob¬ 
tained  by  first  initializing  the  domain  to  the  freestream.  Next,  the  flowfield  was 
iterated  until  the  RR  solution  converged.  Compression  angles  greater  than  or  equal 
to  28°  would  not  resolve  a  RR,  converging  instead  to  a  MR. 

Mach  reflection  solutions  for  9  <  28°  were  obtained  from  a  28°  MR  solution. 
The  28°  MR  solution  was  used  as  the  initial  condition  for  the  27°  MR  solution.  Once 
the  27°  MR  had  converged,  the  resulting  MR  solution  was  used  as  the  starting  point 
for  obtaining  the  26°  MR.  Each  lower  9  used  the  same  procedure  initializing  the 
domain  to  the  converged  MR  solution  of  the  next  higher  9. 

Accelerating  and  decelerating  Lorentz  forces  are  applied  to  the  converged  MR 

and  RR  solutions,  respectively,  as  shown  in  Figures  3.5  and  3.6.  The  applied  non- 
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dimensional  magnetic  field  is  given  a  magnitude  of  1  since  the  magnetic  flux  density, 
B ,  is  non-dimensionalized  by  B0.  The  magnetic  field  for  the  accelerating  Lorentz 
force  is  shown  in  Fig.  3.6  with  B  in  the  —  y  direction.  For  the  decelerating  force,  the 
magnetic  field  is  aligned  with  +y  direction  as  shown  in  Fig.  3.5.  Both  accelerating 
and  decelerating  forces  require  the  electric  field  applied  in  the  positive  ^-direction 
with  E0  nominally  set  to  1.2. 


Figure  3.5  Schematic  of  the  decelerating  Lorentz  force  on  the  25°  x  25°  Regular 
Reflection 

The  Lorentz  force  is  applied  in  a  select  area  defined  by  the  conductivity  tensor 
as  shown  in  Fig.  3.7.  The  electrical  conductivity  tensor,  a,  is  calculated  as  products 
of  modified  Gaussians  of  the  form [4]: 


a  i  =  Ae~a{s-So)n  (3.25) 

where  A  is  set  to  unity,  n  is  set  to  6,  a  is  the  half-width,  Sx  or  5y,  s0  is  the  center  of 
conductivity,  x0  or  y0,  and  s  —  s0  represents  the  distance  from  center  of  conductivity, 
illustrated  in  Fig.  3.7.  Beyond  the  region  defined  by  xq  ±  Sx  and  yo  ±  5y  the  conduc- 
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Figure  3.6  Schematic  of  the  accelerating  Lorentz  force  on  the  26°  x  26°  Mach 
Reflection 

tivity  is  zero.  This  specification  drives  the  electrical  current,  j ,  to  zero  by  Eqn.  2.31, 
subsequently  setting  all  electromagnetic  effects  to  zero. 

The  various  conductivity  specifications  are  shown  in  Table  3.1.  All  computa¬ 
tions  centered  the  applied  force  on  the  centerline  of  the  inlet,  yo=  0.0.  Figure  3.8 
shows  the  six  cases  with  xq  =  0.8175,  which  places  the  center  of  conductivity  co¬ 
incident  with  the  26°  Mach  reflection.  The  conductivity  was  varied  lengthwise  and 
across  the  width  including  Case  2.4.4  with  constant  a. 

Figure  3.9  shows  the  two  cases  with  xo  =  0.4.  At  this  location,  the  applied 
MGD  force  remains  entirely  in  front  of  the  26°  MR.  This  accelerates  the  flow  prior 
to  the  shocks,  with  either  a  small  region  of  force  in  the  center,  Case  1.1.1,  or  a  force 
that  gradually  tapers  away  from  the  center,  Case  1.2.1. 

Figure  3.10  shows  the  two  cases  with  x0  =  1.6.  The  conductivity  is  specified 
to  remain  entirely  behind  a  Mach  or  regular  reflection.  The  flow  behind  the  shock 
formation  is  accelerated  by  the  small  region  of  force  in  the  center,  Case  3.1.1,  or  the 
force  that  gradually  tapers  away  from  the  center,  Case  3.2.1. 
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Figure  3.7 


Area  of  applied  Lorentz  force  defined  by  Equation  3.25 


Table  3.1  Specified  conductivity,  a,  parameters 


Conductivity  Settings 

Case  # 

x0 

yo 

sx 

Sy 

Front  Location,  Center  Focus 

Case  1.1.1 

0.4 

0.0 

0.4 

0.3 

Front  Location,  Gradient  Width 

Case  1.2.1 

0.4 

0.0 

0.4 

0.7 

MR  Location,  Center  Focus 

Case  2.1.1 

0.8175 

0.0 

0.4 

0.3 

MR  Location,  Gradient  Width 

Case  2.2.1 

0.8175 

0.0 

0.4 

0.7 

MR  Location,  Gradient  Full 

Case  2.2.2 

0.8175 

0.0 

1.5 

0.7 

MR  Location,  Full  Width 

Case  2.3.1 

0.8175 

0.0 

50 

50 

MR  Location,  Full  Width,  Half  Length 

Case  2.3.3 

0.8175 

0.0 

0.2 

50 

MR  Location,  Constant  a 

Case  2.4.4 

0.8175 

0.0 

50 

50 

Rear  Location,  Center  Focus 

Case  3.1.1 

1.6 

0.0 

0.4 

0.3 

Rear  Location,  Gradient  Width 

Case  3.2.1 

1.6 

0.0 

0.4 

0.7 
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MR  Location,  Full  Width,  Half  Length 


MR  Location,  Constant  a 


Figure  3.8  Conductivity  centered  at  the  26°  MR,  Cases  2.1.1,  2.2.1  ,  2.2.2,  2.3.1, 
2.3.3,  2.4.4 
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Figure  3.9  Conductivity  centered  in  front  of  the  26°  MR,  Cases  1.1.1,  1.2.1 


Figure  3.10  Conductivity  centered  behind  the  26°  MR  and  the  moved  MR,  Cases 
3.1.1,  3.2.1 
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4-  Results 

The  regular  and  Mach  reflection  solutions  for  compression  angles  of  25°  and  26°  in 
the  dual-solution  domain  were  used  in  this  numerical  study  to  investigate  the  effect 
magnetogasdynamics  has  on  these  shock  structures.  Figures  1.5  and  1.6  show  the  25° 
RR  and  MR  numerical  solutions,  respectively.  Figure  4.1  shows  the  26°  dual-solution 
domain  with  the  RR  presented  on  bottom  and  the  MR  on  top.  For  this  numerical 
study,  the  25°  RR  was  subjected  to  a  decelerating  Lorentz  force  in  order  to  determine 
the  sensitivity  of  transitioning  a  RR  to  an  undesired  MR.  The  26°  MR  was  subjected 
to  an  accelerating  Lorentz  force  to  determine  its  capability  of  transitioning  a  MR  to 
a  desired  RR  and  suppressing  inlet  unstart. 


Figure  4.1  Pressure  contours  from  a  regular  and  Mach  reflection  for  Mach  4.96, 
2-D  solutions,  26°  x  26°  (MR-right,  RR-left) 

The  MGD  code  uses  temperature,  Mach  number,  and  Reynolds  number  to 
calculate  the  flow  variables  for  each  solution.  The  current  research  used  freestream 
flow  parameters  of  Re  =  2  x  106  (based  on  w  —  1  m),  M0  =  4.96,  T0  =  72.54  K,  p0  = 
1.1614x  10-2  kg/m 3,  po  =  241.776  Pa,  Vo  =  846.76  m/s,  and  t0  =  1.118x  10~3  s  were 
employed  in  this  research.  The  corresponding  total  pressure  and  total  temperature 
were  1.22  x  105  Pa  and  429.43  K ,  respectively.  These  parameters  were  chosen  to 
match  the  experimental  work  of  Ivanov,  et  al.  and  the  numerical  work  of  Schmisseur 
and  Gaitonde[17,  20,  29,  30]. 
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The  magnetogasdynamic  effects  were  applied  to  the  flowfield  using  a  magnetic 
flux  density  set  to  a  nominal  value  of  B0  =  2.0  T  for  Q  =  1.  The  resulting  maximum 
magnitude  of  the  conductivity,  cr0,  was  2.459  mho/m  based  on  Eqn.  2.42.  The 
magnitude  of  the  applied  electric  field,  Eq,  is  2032.23  v/m  from  Eqn.  2.43.  These 
values  resulted  in  a  magnetic  Reynolds  number,  Rm  =  0.0026  (Eqn.  2.37),  much 
less  than  1  as  required  for  the  source  term  formulation  of  the  MGDE.  When  the 
interaction  parameter  was  increased,  the  magnitude  of  the  conductivity  was  kept 
constant  and  the  magnetic  field  was  increased  as  required.  For  this  research,  Q  = 
10.0  was  the  largest  interaction  parameter  investigated,  which  equates  to  a  required 
magnetic  flux  density  magnitude  of  6.325  T  and  a  maximum  applied  electric  field  of 
6426.49  v/m. 

4.1  Hysteresis 

Schmisseur  and  Gaitonde’s  hysteresis  for  g/w  =  0.34  (Fig.  1.10)  [30]  is  dupli¬ 
cated  in  Fig.  4.2.  The  shock  structure  remains  in  a  RR  configuration  until  9  >  27°. 
After  transitioning  to  a  MR  configuration,  the  Mach  reflection  remains  until  9  <  22°. 
This  hysteresis  matches  Schmisseur  and  Gaitonde’s  numerical  results [30],  the  initial 
point  for  the  magnetogasdynamic  study  in  this  research.  This  numerical  work  also 
corresponds  to  the  experimental  work  of  Ivanov,  et  al.  [17,  20] 

4-2  Grid  Convergence 

A  grid  convergence  study  was  performed  to  show  that  the  grids  used  in  this 
research  have  a  high  enough  resolution  to  accurately  capture  the  fluid  interactions. 
The  26°  Mach  reflection  solution  was  used  to  compare  the  coarse,  fine,  and  ultra-fine 
grids  shown  in  Fig.  3.2.  Figure  4.3  shows  the  coarse  mesh  converged  MR  on  bottom 
and  the  fine  mesh  converged  MR  on  top.  The  increased  cell  density  resolves  a  normal 
shock  on  the  fine  mesh  that  is  40.3%  more  compact  than  its  coarse  mesh  counterpart. 
As  illustrated  in  the  enlarged  view,  the  normal  shocks  fall  within  one  coarse  cell  or 
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Increasing  Wedge  Angle,  6 


Figure  4.2  Two-dimensional  dual-solution  domain,  g/w  =  0.34,  =  4.96 


two  fine  cells  of  each  other.  The  center  of  the  coarse  normal  shock  is  1.72  coarse 
cells  away  from  the  corresponding  location  on  the  fine  normal  shock.  This  distance 
is  a  2.44%  movement.  The  Mach  stem  height  increases  by  2.54%  from  the  coarse  to 
the  fine  solution.  The  high  correlation  between  these  solutions  validates  the  use  of 
the  coarse  mesh  for  preliminary  work  and  trend  determination. 

The  slight  difference  between  the  coarse  and  fine  mesh  fluid  interaction  resolu¬ 
tion  is  enough  to  change  the  flow  behind  a  regular  reflection.  Figure  4.4  shows  the 
converged  RR  for  both  grids.  On  bottom,  the  coarse  mesh  resolves  the  primary  and 
secondary  oblique  shocks,  the  expansion  fan,  the  shock-expansion  fan  interactions, 
and  the  oblique  shock  reflections  continuing  downstream  in  the  inlet,  as  expected. 
However,  on  top  the  solution  resolved  on  the  fine  mesh  includes  an  unexpected  fea- 
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Figure  4.3  MR  Pressure  contours  for  coarse,  107  x  50  on  bottom,  and  fine,  213  x  99 
on  top,  grids,  26°  x  26° 

tnre.  Instead  of  traversing  the  expansion  fan  and  reflecting  off  the  wall,  the  reflected 
oblique  shock  forms  a  small  Mach  reflection  near  the  wall. 

Figure  4.5  shows  u- velocity  contours  for  this  26°  RR.  The  u- velocity  remains 
positive  everywhere  on  the  coarse  mesh,  but  a  separation  region  behind  the  secondary 
MR  is  present  in  the  flowfield  on  the  fine  mesh.  The  enlarged  view  in  Fig.  4.5  shows 
this  separation  region.  This  flow  phenomena  is  due  to  the  induced  curvature  of  the 
reflected  shock  caused  by  the  shock-expansion  fan  interaction.  The  velocity  near 
the  wall  is  accurately  resolved  by  the  fine  mesh  due  to  the  greater  cell  density  in 
the  region  near  the  wall.  As  the  solution  converges,  the  angle  of  the  flow  is  turned 
perpendicular  to  the  wall.  The  wall  slip  BC  sets  the  velocity  tangent  to  the  wall. 
In  this  case,  as  the  flow  is  turned  perpendicular,  the  velocity  on  the  wall  decreases. 
When  the  velocity  becomes  subsonic,  a  shock  develops.  Once  the  shock  is  formed,  a 
larger  region  of  flow  is  reduced  to  subsonic  velocities.  In  order  to  support  this  increase 
in  subsonic  flow,  the  shock  moves  forward  causing  a  separation  to  occur.  The  pressure 
change  across  the  oblique  shocks  grows  too  large,  and  a  Mach  reflection  develops. 
This  phenomena  does  not  develop  on  the  coarse  mesh  due  to  the  lower  resolution. 
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Secondary  Mach  Reflection 


Figure  4.4  RR  pressure  contours  for  coarse  and  fine  grids  on  bottom  and  top, 
respectively,  26°  x  26° 

In  addition,  this  secondary  Mach  reflection  only  appears  in  regular  reflections  where 

e  >  26°. 

Figure  4.6  shows  the  fine  mesh  converged  MR  on  bottom  and  the  ultra-fine 
mesh  converged  MR  on  top.  The  ultra- fine  normal  shock  is  51.2%  more  compact  than 
the  normal  shock  resolved  on  the  fine  mesh.  The  center  of  the  fine  normal  shock  is 
1.2  fine  cells  away  from  the  center  of  the  ultra- fine  normal  shock,  corresponding  to  a 
movement  of  1.07%.  The  enlarged  view  depicts  normal  shocks  that  virtually  overlap. 
Figure  4.6  shows  good  agreement  between  the  fine  and  ultra-fine  mesh  solutions  with 
the  Mach  stem  height  increasing  by  0.87%  from  the  fine  to  ultra-fine  solution.  Since 
the  ultra-fine  mesh  required  approximately  five  times  longer  to  converge  than  the 
fine  mesh,  the  fine  mesh  was  chosen  as  the  primary  mesh  for  this  numerical  study.  A 
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Magnified  View 
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Figure  4.5  RR  u-velocity  contours  for  coarse  and  fine  grids,  on  bottom  and  top, 
respectively,  26°  x  26°  (Fine-Magnified) 


fine  mesh  was  used  for  all  solutions  presented  through  the  rest  of  this  section  unless 
otherwise  noted. 


4-3  Decelerating  Lorentz  Force 

In  order  to  determine  the  regular  reflection  flowfield’s  sensitivity  to  transition  to 
Mach  reflection,  the  25°  RR  flowfield  shown  in  Fig.  1.5  was  subjected  to  a  decelerating 
Lorentz  force  as  defined  in  Fig.  3.5.  The  MGD  force  was  specified  by  the  conductivity 
pattern  of  case  2.1.1(Fig.  3.8),  a  small  bread  loaf  pattern  focused  over  the  26°  MR. 
To  accomplish  this  deceleration,  the  magnetic  held  was  applied  in  the  +y  direction 
and  the  interaction  parameter  was  set  to  1.0.  Figure  4.7  shows  the  solutions  as  the 
decelerating  force  was  applied  for  a  non-dimensional  time,  r,  of  0.02,  0.05,  0.1,  0.2, 
and  0.3  which  correspond  to  2.36  x  10“5  s,  5.9  x  10~5  s,  1.18  x  10_4  s,  2.36  x  10_4  s, 
and  3.54  x  10”4  s,  respectively.  The  time  step,  At,  was  1.18  x  10~6  s,  At  =  0.001. 
Each  inlet  represents  a  snapshot  of  the  flowfield  at  the  specified  r  with  the  force 
applied.  In  Fig.  4.7,  the  bottom  half  of  each  inlet  is  the  25°  RR  and  the  top  half  shows 
the  effect  of  the  force  on  the  RR.  Force  applications  less  than  r  =  0.02  have  no  visible 
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Figure  4.6  MR  Pressure  contours  for  fine,  213  x  99  on  top,  and  ultra-fine,  425  x  197 
on  bottom,  grids,  26°  x  26° 

effect  on  the  the  initial  and  reflected  oblique  shocks.  A  decelerating  Lorentz  force 
applied  for  r  =  0.05  moves  the  reflection  point  of  oblique  shocks  upstream,  beginning 
the  transition  to  a  normal  shock  in  the  middle  of  the  inlet.  This  movement  is  due 
to  the  decrease  in  the  local  Mach  number  to  3.12.  The  lower  local  Mach  number 
increases  the  oblique  shock  angle  f3.  As  the  time  of  force  application  increases,  the 
pressure  contours  move  upstream  and  form  a  normal  shock.  The  pressure  contours 
of  the  r  =  0.3  solution  are  aligning  at  the  location  of  the  25°  MR  as  shown  in 
Fig.  1.6.  This  decelerating  Lorentz  force  destabilizes  the  regular  reflection  enough 
in  2.36  x  10  5  s  to  transition  the  desirable  RR  shock  system  to  the  inefficient  Mach 
reflection. 

Figure  4.8  compares  two  solutions  with  the  decelerating  Lorentz  force  applied 
for  r  =  0.02.  The  top  solution  is  200  iterations  of  the  decelerating  force  with 
a  At  =  0.0001  and  the  bottom  solution  is  the  same  force  after  20  iterations  of 
At  =  0.001.  The  excellent  agreement  between  the  two  solutions  shows  that  this 
MGD  application  is  independent  of  time  step. 
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t  =  0.02 


t  =  0.05 


Figure  4.7  Comparison  of  the  RR  movement  due  to  a  decelerating  Lorentz  force 
with  increasing  r  (RR-bottom,  MGD-top) 

When  the  results  of  the  decelerating  force  were  subjected  to  closer  examination, 
the  sensitivity  of  the  regular  reflection  to  transition  was  easily  observed.  The  pressure 
contours  of  the  RR  have  moved  upstream  when  r  >  0.05.  This  movement  of  the 
shock  system  signifies  the  beginning  of  the  transition  to  a  MR  as  shown  in  Fig.  4.7. 
Decelerating  Lorentz  force  applications  less  than  2.36  x  10_5  s  did  not  destabilize  the 
RR,  and  the  regular  reflection  shock  configuration  remained  even  after  the  force  was 
removed.  For  force  applications  of  r  >  0.02,  but  less  than  or  equal  to  0.7,  the  shock 
system  transitions  to  Mach  reflection  when  the  force  is  removed.  These  MRs  formed 
by  the  decelerating  Lorentz  force  compare  favorably  with  the  25°  MR.  Figure  4.9 
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Figure  4.8  Decelerating  Lorentz  force  Time  Step  comparison  pressure  contours, 
25°  x  25° 

shows  the  25°  converged  MR  on  bottom  and  the  MR  formed  by  a  decelerating  Lorentz 
force  applied  for  r  =  0.1  on  top.  The  regular  reflection  is  extremely  sensitive  to  the 
disturbance  created  by  this  decelerating  Lorentz  force.  The  inlet  sustains  higher 
total  pressure  losses  and  a  lower  compression  ratio  due  to  this  RR  transition  to  MR, 
and  a  higher  probability  of  inlet  unstart. 

The  inlet  unstarts  when  the  decelerating  Lorentz  force  is  applied  to  the  regular 
reflection  for  longer  than  8.267  x  10-4  s,  r  =  0.7.  The  rate  of  this  inlet  unstart  is 
dependent  upon  how  long  the  force  is  applied.  Figure  4.10  demonstrates  the  effects 
of  the  decelerating  force  when  the  application  period  was  increased  beyond  those 
shown  in  Fig.  4.8.  The  25°  MR  is  shown  on  the  bottom  of  each  inlet  and  the 
decelerating  Lorentz  force  solution  for  each  r  is  shown  on  top.  When  the  force  is 
applied  for  r  >  0.3,  the  shock  structure  is  moved  upstream  beyond  the  25°  MR. 
A  decelerating  Lorentz  force  applied  longer  than  8.267  x  10~4  s  created  a  normal 
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Figure  4.9  MR  formed  by  a  decelerating  Lorentz  Force  (25°  MR-bottom,  MGD 
MR-top) 

shock  that  propagates  to  the  entrance  of  the  inlet.  This  shock  system  was  expelled 
even  if  the  force  was  removed.  The  inlet  unstart  dramatically  reduces  the  thrust  of 
the  SCRAMjet  and  increases  the  drag.  Figure  4.11  illustrates  inlet  unstart  due  to  a 
decelerating  Lorentz  force.  The  force  was  applied  for  r  =  1.5,  1.771  x  ICR3  s,  and 
then  removed.  The  inlet  unstarts  as  the  flowfield  converges. 

The  regular  reflection  is  extremely  sensitive  to  any  upstream  disturbance,  tran¬ 
sitioning  to  a  Mach  reflection  almost  immediately.  Transient  application  of  the  de¬ 
celerating  Lorentz  force  for  microseconds  transitioned  the  RR  to  the  inefficient  MR 
shock  configuration.  The  higher  total  pressure  losses  and  reduced  compression  ratio 
of  this  MR  configuration  will  reduce  the  thrust  of  a  SCRAMjet.  Applications  of  this 
force  for  milliseconds  resulted  in  inlet  unstart. 
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Increasing  decelerating 
force  application  time 


Figure  4.10  Comparison  of  the  RR  movement  due  to  a  decelerating  Lorentz  force 
with  increasing  r  (25°  MR-bottom,  MGD-top) 


4-4  Accelerating  Lorentz  Force 

This  numerical  study  investigates  the  use  of  the  interaction  between  electro¬ 
magnetic  fields  and  a  conductive  flowfield  to  reduce  the  inefficiencies  in  a  SCRAMjet 
inlet  created  by  the  transition  of  a  regular  reflection  to  a  Mach  reflection.  The  26° 
MR  was  subjected  to  an  accelerating  Lorentz  force  with  the  goal  of  reducing  these 
inefficiencies  and  possibly  transitioning  the  MR  to  a  RR.  The  accelerating  Lorentz 
force  was  created  by  applying  a  magnetic  field  in  the  —y  direction  and  an  electric 
field  in  the  +z  direction  as  shown  in  Fig.  3.6.  This  MGD  force  was  applied  in  the 
region  specified  by  the  conductivity  (Eqn.  3.25).  The  region  of  conductivity  was 
varied  through  all  of  the  cases  presented  in  Table  3.1  and  shown  in  Figures  3.8,  3.9, 
and  3.10.  These  cases  were  numerically  implemented  on  the  coarse  mesh  to  decrease 
the  computational  time  required  in  determining  the  optimum  conductivity  case. 
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Figure  4.11  Inlet  unstart  due  to  a  decelerating  Lorentz  Force 

The  area  where  the  accelerating  Lorentz  force  acts  is  determined  entirely  by 
the  conductivity  pattern.  A  pattern  completely  in  front  of  the  Mach  reflection  only 
accelerates  the  fluid  before  it  passes  through  the  normal  shock.  The  patterns  entirely 
behind  the  Mach  reflection  accelerate  the  fluid  downstream  of  the  normal  shock,  and 
the  fluid  upstream  of  the  shock  system  is  not  influenced  by  this  acceleration.  In  both 
cases,  one  would  expect  the  limited  application  of  the  accelerating  force  to  have  very 
little  effect  on  the  shock  system.  Centering  the  accelerating  Lorentz  force  over  the 
MR  accelerates  the  fluid  upstream  and  downstream  of  the  normal  shock.  The  MR 
location  should  push  the  MR  downstream  further  and  reduce  the  Mach  stem  more 
than  either  of  the  other  two  application  locations.  Applying  the  force  over  greater 
areas  of  the  flowfield  in  both  length  and  width  increases  the  energy  addition  to 
the  flow.  This  greater  energy  addition  should  result  in  greater  acceleration  and 
greater  MR  movement.  The  optimum  conductivity  pattern  should  have  the  greatest 
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downstream  movement  of  the  normal  shock  without  inducing  flowfield  distortion, 
such  as  the  secondary  MRs. 

The  optimum  a  pattern  was  implemented  and  the  interaction  parameter  was 
varied  from  1.0  to  10.0  to  determine  its  effect  on  the  MR  flowfield.  The  increases  in  Q 
represent  increases  in  the  applied  magnetic  field.  The  source  term  on  the  continuity 
of  momentum  equation  is  Q(j  x  D)  and  the  magnitude  of  the  applied  magnetic  flux 
is  defined: 


B0  = 


/ 


Qp()V{) 


(4.1) 


The  current,  j.  increases  linearly  with  the  magnetic  flux.  The  combination  of  these 
effects  specifies  that  the  magnitude  of  the  Lorentz  force  goes  with  Q2.  Increases 
in  the  interaction  parameter  should  result  in  a  more  effective  force,  magnifying  the 
changes  created  in  the  flowfield. 


The  load  factor  was  varied  independent  of  the  interaction  parameter  to  de¬ 
termine  its  effect  on  the  MR.  /C  was  increased  from  the  nominal  values  of  1.2,  to 
1.8  and  then  to  2.4,  increases  of  50%  and  100%,  respectively.  The  increases  in  /C 
represent  increases  in  the  applied  electric  field.  The  force  should  increase  linearly 
with  the  load  factor  because  the  current  increases  linearly  with  the  electric  field. 
This  increase  is  rapidly  offset  by  the  square  relationship  between  the  heating  and  /C. 
The  increases  in  the  load  factor  should  actually  be  detrimental  to  the  flowfield.  The 
greater  Lorentz  force  will  not  be  able  to  offset  the  greater  heating  and  the  MR  will 
move  upstream,  or  the  inlet  could  possibly  unstart. 


The  mass-averaged  static  pressure  increases  by  a  factor  of  22.18  through  the 
shock  structure  of  a  26°  regular  reflection  at  =  4.96.  The  Mach  reflection  for 
identical  conditions  only  compresses  the  flow  by  a  factor  of  19.68,  11%  less  than  the 
RR  compression  ratio.  The  mass-averaged  total  pressure  ratios  across  26°  RR  and 
MRs  are  0.283  and  0.249,  respectively.  The  transition  from  RR  to  MR  also  moves 


4-13 


the  shock  system  upstream  23%  to  a  less  stable  location  in  front  of  the  throat.  Use 
of  an  accelerating  Lorentz  force  to  trigger  the  transition  to  a  RR  is  desirable  because 
of  this  12%  loss  in  total  pressure  recovery  and  11%  reduction  in  compression.  This 
force  would  reduce  the  total  pressure  losses  associated  with  the  normal  shock  of  a 
MR  and  provide  better  flow  compression  with  the  transitioned  RR.  Additionally, 
this  force  has  the  capability  to  prevent  potential  engine  inlet  unstart  due  to  the 
instability  of  the  MR. 

4-4- 1  Specification  of  Conductivity.  The  coarse  mesh  is  used  to  demon¬ 
strate  the  overall  effects  of  the  ten  conductivity  patterns  used  in  this  research  (Table  3.1). 
While  the  results  are  less  accurate  than  those  on  the  fine  mesh,  the  general  trends 
are  duplicated.  Fewer  nodes  decrease  the  computational  time  required  for  conver¬ 
gence.  For  this  numerical  study,  the  conductivity  pattern  that  provides  the  largest 
reduction  in  the  Mach  reflection  without  overly  distorting  the  flowheld  is  the  most 
desirable  pattern.  The  feasibility  of  creating  this  held  is  not  examined. 

Figure  4.12  illustrates  the  effect  that  changing  the  conductivity  specification 
has  on  the  flowheld.  The  26°  MR  is  shown  on  the  bottom  of  each  inlet  and  the 
howheld  with  the  Q  =  3.0  accelerating  force  is  presented  on  top.  The  focused 
conductivity,  case  2.1.1,  has  a  very  large  effect  on  the  MR  pushing  the  normal  shock 
downstream  to  a  position  that  is  34.57%  of  the  difference  between  the  26°  Mach 
and  regular  reflection  locations.  The  Mach  stem  height  is  reduced  also.  Application 
of  this  conductivity  pattern  also  destroys  the  reflection  pattern  further  downstream 
developing  a  new  Mach  reflection  near  the  wall.  These  flowheld  effects  are  the  result 
of  the  how  directly  in  the  center  of  the  inlet  accelerating  with  little  or  no  change  in 
the  how  outside  the  a  concentration,  creating  effects  much  like  a  jet. 

In  case  2.2.1  the  conductivity  is  focused  directly  over  the  MR  and  is  gradually 
reduced  near  the  walls.  Figure  4.12  illustrates  that  this  a  pattern  results  in  a  MR 
moved  further  downstream  than  the  case  2.1.1,  reducing  the  distance  between  the 
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Figure  4.12  Contours  of  pressure  with  Q  =  3.0,  conductivity  centered  at  the  26° 
MR,  Cases  2.1.1,  2.2.1  ,  2.2.2,  2.3.1,  2.3.3,  2.4.4,  Coarse  Mesh  (MGD- 
top,  MR-bottom) 
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MR  and  RR  by  84%.  The  greater  movement  is  accomplished  without  the  flowfield 
degeneration  of  case  2.1.1.  The  flowfield  when  o  spans  the  entire  width  of  the  inlet, 
2.3.1,  is  very  similar  to  2.2.1.  The  normal  shock  is  moved  almost  as  far  downstream 
to  a  location  that  is  16.87%  of  the  original.  The  reduced  length  specification  in  case 
2.3.3  greatly  decreases  the  impact  of  the  accelerating  force,  only  moving  the  normal 
shock  downstream  30%.  These  cases  with  greater  width  specification  do  not  create 
the  jet  effect  of  case  2.1.1  because  the  flow  all  the  way  across  the  inlet  is  accelerated. 
The  reduced  length  of  the  2.3.3  pattern  accelerates  the  flow  for  a  shorter  period  and 
correspondingly  has  less  effect. 

Cases  2.2.2  and  2.4.4,  the  gradient  full  and  constant  a  specifications  respec¬ 
tively,  have  the  greatest  impact  on  the  Mach  reflection  as  shown  in  Fig.  4.12.  The 
MR  is  moved  to  a  location  0.83%  downstream  of  the  26°  RR,  a  100.83%  movement 
for  case  2.2.2.  In  this  case  the  conductivity  is  a  constant  1.0  for  most  of  the  flowfield 
and  gradually  tapers  down  at  the  entrance  and  exit  of  the  inlet.  A  constant  a  —  1.0 
throughout  the  flowfield,  case  2.4.4,  moves  the  MR  downstream  to  a  location  95.8% 
of  it  original  location.  The  accelerating  force  in  both  cases  induces  none  of  the  added 
flowfield  distortion  present  in  the  center  focused  pattern,  case  2.1.1.  The  flowfield 
in  both  cases  is  uniformly  accelerated  in  width  and  length,  accelerating  the  flowfield 
more  than  the  other  cases. 

The  effects  of  the  accelerating  Lorentz  force  are  very  similar  for  cr  cases  2.2.2 
and  2.4.4.  Figure  4.13  shows  the  gradient  full  pattern  on  bottom  and  the  constant 
a  on  top.  The  gradient  full  conductivity  pattern  moves  the  normal  shock  further 
downstream  than  the  constant  a.  This  movement  corresponds  to  a  smaller  Mach 
stem.  While  case  2.2.2  and  2.4.4  are  similar,  Fig.  4.14  illustrates  that  the  gradient 
specification  is  the  more  optimum  of  the  two.  Case  2.2.2  is  on  the  bottom  of  the 
inlet  and  case  2.4.4  is  on  top.  The  constant  a  pattern  has  a  greater  maximum 
temperature,  and  greater  temperature  rise  overall,  than  the  gradient  full  pattern. 
The  gradient  pattern  adds  less  energy  to  the  flowfield.  As  soon  as  the  flow  enters 
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the  inlet,  the  constant  a  pattern  adds  the  maximum  energy  to  the  flow.  This  high 
input  shows  up  in  the  temperature  rise. 


Figure  4.13  Contours  of  pressure  with  Q  =  3.0,  Case  2.2.2  on  bottom  and  case 
2.4.4  on  top  on  the  coarse  Mesh 

Figure  4.15  shows  the  effects  of  the  accelerating  force  when  the  force  is  applied 
in  front  of  the  MR.  Cases  1.1.1  and  1.2.1  are  shown  at  the  top  of  the  inlets  and  the 
26°  MR  is  on  the  bottom.  Both  conductivity  patterns  effect  the  flowfield,  but  the 
effects  of  the  accelerating  Lorentz  force  are  less  than  the  MR  centered  patterns  when 
the  force  is  applied  in  front  of  the  Mach  reflection.  Case  1.1.1  moves  the  normal 
shock  downstream  22.45%  and  case  1.2.1  moves  the  normal  shock  to  a  location  that 
is  19.93%  less  than  the  original  value.  The  Mach  stem  height  is  reduce  accordingly 
for  both  cases.  The  focused  pattern,  case  1.1.1,  distorts  the  flowfield  slightly  more, 
as  expected,  due  to  the  center  focus  of  the  accelerating  force  like  case  2.1.1. 

Figure  4.16  shows  the  effects  of  the  accelerating  force  when  it  is  applied  behind 
the  MR.  The  accelerating  MGD  force  is  depicted  at  the  top  of  the  inlets  and  the  26° 
MR  is  on  the  bottom.  The  rear  location  has  very  little  effect  on  the  shock  system  for 
either  pattern,  moving  the  MR  back  12.34%  in  both  cases.  These  minimal  changes 

are  due  to  application  of  the  accelerating  force  behind  the  normal  shock  in  both 
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Figure  4.14  Contours  of  temperature  with  Q  =  3.0,  Case  2.2.2  on  bottom  and  case 
2.4.4  on  top  on  the  coarse  Mesh 


Figure  4.15  Contours  of  pressure  with  Q  =  3.0,  conductivity  centered  in  front  of 
the  26°  MR,  Cases  1.1.1,  1.2.1,  Coarse  Mesh  (MGD-top,  MR-bottom) 
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conductivity  patterns.  The  supersonic  flow  upstream  of  the  shock  system  can  not 
respond  to  the  downstream  application. 


Figure  4.16  Contours  of  pressure  with  Q  =  3.0,  conductivity  centered  behind  the 
26°  MR,  Cases  3.1.1,  3.2.1,  Coarse  Mesh  (MGD-top,  MR-bottom) 

The  effects  of  the  conductivity  pattern  are  independent  of  the  interaction  pa¬ 
rameter.  The  gradient  full  conductivity  pattern,  case  2.2.2,  has  a  greater  effect  on 
the  Mach  reflection  and  less  flowfield  distortion  than  the  focused  conductivity,  case 
2.1.1,  regardless  of  the  value  of  Q.  Figure  4.17  compares  these  two  conductivity  pat¬ 
terns  at  interaction  parameters  of  1.0,  3.0,  and  5.0.  All  inlets  have  the  gradient  full 
pattern  on  top  and  the  focused  pattern  on  bottom.  For  the  equivalent  interaction 
parameter,  the  Mach  reflection  is  moved  much  further  downstream  by  the  gradient 
full  pattern.  The  temperature  contours  illustrate  the  severe  problems  created  by  the 
center  focused  a  pattern.  Case  2.1.1  with  Q  =  3.0  creates  non-dimensional  temper¬ 
atures  in  excess  of  1850.0.  In  order  to  show  any  of  the  temperature  contours  for 
case  2.2.2  with  Q  =  3.0,  the  temperature  contours  for  the  whole  Q  =  3.0  inlet  are 
logarithmic.  With  the  linear  scale  used  on  the  other  inlets  only  one  contour  would 
appear  for  case  2.2.2  if  30  temperature  contours  were  used  because  its  maximum 
non-dimensional  temperature  is  68.67.  The  maximum  temperature  for  case  2.1.1 
when  Q  =  3.0  is  equivalent  to  a  temperature  in  excess  of  135,000  K ,  well  outside  the 
bounds  of  the  ideal  gas  law.  This  temperature  is  beyond  the  bounds  of  the  models 
used  in  this  research,  but  the  general  trend  is  correct  and  the  heating  encountered 

using  a  focused  conductivity  is  excessive.  The  largest  amount  of  energy  added  to  the 
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flow  occurs  along  the  centerline.  By  the  end  of  the  region  of  conductivity,  the  Ae 
along  the  center  of  the  inlet  is  at  the  maximum  value.  The  Q  =  3.0  focused  a  MR 
is  located  at  this  point  of  maximum  Ae.  The  combination  of  the  increased  heating 
due  to  the  added  energy  and  the  temperature  rise  across  the  normal  shock  results  in 
the  maximum  temperature  encountered  in  this  research.  Heating  excesses  such  as 
this  are  outside  the  models  used  in  the  code  employed  in  this  research. 


Figure  4.17  Comparison  of  the  MR  movement  due  to  an  accelerating  Lorentz  force 
with  increasing  interaction  parameter  values,  a  Case  2.1.1  on  bottom 
and  case  2.2.2  on  top,  Fine  Mesh 


When  Q  =  5.0,  the  temperatures  for  case  2.1.1  are  well  within  the  range  of  the 
ideal  gas  law  with  a  maximum  value  of  roughly  2000  K.  While  these  temperatures 
are  outside  of  the  constant  7  bounds,  the  area  where  the  temperature  is  greater  than 
800  K  is  very  small,  and  the  presented  trends  accurate.  The  MGD  code  employed 
for  this  numerical  study  is  valid  for  most  of  the  flowfield,  but  future  work  might 
include  the  capability  for  a  non-constant  ratio  of  specific  heats. 

The  total  pressure  losses  are  greater  for  case  2.1.1  than  case  2.2.2  as  shown  in 
Fig.  4.18.  The  gradient  full  pattern  is  again  presented  on  the  top  of  the  inlet  and  the 
focused  pattern  is  on  the  bottom.  The  average  Mach  number  is  greater  when  the 
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gradient  full  a  pattern  is  used.  The  separation  region  and  secondary  Mach  reflection 
present  in  case  2.1.1  are  not  present  in  case  2.2.2  on  the  coarse  or  fine  mesh. 


Figure  4.18  Q  =  5.0  parallel  Lorentz  force  Mach  and  total  pressure  contours,  o 
Case  2.1.1  on  bottom  and  case  2.2.2  on  top,  Fine  Mesh 


The  gradient  full  conductivity  pattern  is  the  most  effective  of  the  ten  conductivity 


patterns  studied  in  this  research  and  includes  the  least  flowfield  degradation.  Unless 


otherwise  specified,  case  2.2.2  is  used  for  the  rest  of  this  section. 


4-4-2  Magnetic  Field  Magnitude.  In  order  to  try  to  transition  the  MR  to  a 
RR,  the  converged  26°  MR  solution  was  subjected  to  an  accelerating  Lorentz  force 
with  interaction  parameters  from  1.0  to  10.0  and  the  conductivity  specified  by  case 
2.2.2.  The  load  factor  remained  at  the  nominal  value  of  1.2  specified  by  Eq  =  1.2. 
When  the  accelerating  Lorentz  force  was  acting  parallel  to  the  flow,  the  26°  MR 
moved  downstream  through  the  throat  to  a  new  location  47.55%  of  the  original  MR 
location  when  Q  =  1.0.  Figure  4.19  illustrates  the  changes  in  the  MR  due  to  the 
applied  force.  In  this  figure,  the  26°  RR  solution  is  depicted  on  the  bottom  of  the 
flowfield,  and  the  MR  with  application  of  the  accelerating  Lorentz  force  is  depicted 
on  top.  As  the  interaction  parameter  is  increased  from  1.0  to  10.0,  the  magnitude 
of  the  applied  magnetic  field  is  increased,  and  the  MR  moves  further  downstream. 

While  the  Mach  reflection  is  pushed  beyond  the  location  of  the  regular  reflec¬ 
tion,  the  MR  is  not  fully  transitioned  to  a  RR  reflection  in  Fig.  4.19.  The  shock 
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Increasing  Interaction  Parameter,  Q 

Figure  4.19  Movement  of  the  MR  due  to  a  parallel  Lorentz  force  with  increasing 
interaction  parameter  values  (MGD-top,  RR-bottom) 

systems  of  both  the  Q  =  7.5  and  10.0  solutions  are  pushed  downstream  well  beyond 
the  location  of  the  26°  RR  to  new  locations  that  are  27.74%  and  36.97%  beyond  the 
regular  reflection,  respectively.  The  locations  of  the  new  MRs  are  a  result  of  the  the 
increased  flowheld  Mach  numbers.  By  the  time  the  flow  reaches  the  normal  shock  in 
the  Q  =  10.0  solution,  the  local  Mach  number  is  5.44.  This  increased  Mach  number 
results  in  a  lower  oblique  shock  angle,  (3.  Figure  4.20  depicts  the  Q  =  10.0  MR  on 
top  of  the  each  inlet  with  the  26°  MR  and  RR  on  the  bottom.  The  top  inlets  show 
contours  of  Mach  number,  the  middle,  contours  of  pressure,  and  the  bottom  inlets, 
contours  of  total  pressure.  The  oblique  shock  is  bent  towards  the  wall  due  to  the 
Lorentz  force  acceleration.  The  average  Mach  number  and  total  pressure  at  the  exit 
of  the  inlet  is  higher  in  the  MGD  inlet  than  the  26°  RR  inlet.  This  is  dictated  by  the 
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energy  interaction  and  Lorentz  force  source  terms  in  the  MGDE.  The  accelerating 
Lorentz  force  increases  the  velocity  of  the  flow  in  the  region  specified  by  a.  Total 
pressure  is  based  on  pressure  and  Mach  number,  and  the  Mach  number  increases  as 
the  velocity  increases.  This  inter-relation  would  allow  the  total  pressure  to  increase 
through  the  inlet,  given  a  high  enough  Q.  The  26°  RR  solution  also  experiences 
greater  losses  in  Mach  number  and  total  pressure  due  to  the  secondary  MR. 


Total  Pressure 


Q  =  10.0  vs.  MR 


Mach 


Pressure 


Figure  4.20  Q  =  10.0,  MR,  and  RR  parallel  Lorentz  force  Mach,  pressure,  and 
total  pressure  contours(MGD-top,  MR  &  RR-bottom) 
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Figures  4.21  and  4.22  show  the  regular  and  Mach  reflections  respectively,  on 
the  bottom  of  an  inlet  with  the  Q  =  10  flowfield  on  top.  In  the  enlarged  view, 
the  contours  of  Mach  show  that  the  normal  shock  is  95%  smaller  than  the  26° 
MR,  but  the  higher  local  Mach  number  induced  by  the  accelerating  Lorentz  force 
creates  a  much  stronger  normal  shock.  The  flowfield  behind  the  shock  structure 
when  Q  =  10.0  displays  the  expected  pattern  of  oblique  shock  reflections  unlike  the 
flowfield  behind  the  26°  RR.  The  energy  added  to  the  flowfield  prevents  the  flow 
from  forming  the  secondary  Mach  reflections  present  in  the  RR. 


4.91088 

4.37822 

3.84556 

3.3129 

2.78025 

2.24759 

1.71493 

1.18227 

0.649612 

0.116954 


Figure  4.21  Mach  contours  for  Q  =  10.0  on  top  and  the  26°  RR  on  bottom 


In  order  to  provide  a  measure  of  the  effectiveness  of  an  accelerating  Lorentz 
force,  mass-averaging  was  used.  The  location  where  the  data  sampling  takes  place 
can  greatly  effect  the  results  of  the  mass-averaging.  The  changing  location  of  the 
shock  system  makes  the  sampling  difficult  in  this  research.  The  26°  MR  and  RR  are 
chosen  as  the  baseline  parameters  and  each  solution  is  sampled  behind  the  throat, 
at  a  location  0.12  behind  the  normal  shock  or  reflection. 
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Figure  4.22  Mach  contours  for  Q  =  10.0  on  top  and  the  26°  MR  on  bottom 


Though  transition  does  not  occur,  many  of  the  desired  goals  are  realized  by 
this  application  of  an  accelerating  Lorentz  force.  Table  4.1  lists  the  mass-averaged 
ratios  of  pressure,  temperature,  total  pressure,  and  total  temperature  across  the 
shock  system,  and  includes  the  maximum  temperature  in  the  flowfield,  for  the  RR, 
MR,  and  Q  =  1.0  —  10.0  solutions.  The  RR  and  MR  values  for  pressure  and  tem¬ 
perature  should  not  be  compared  to  the  corresponding  values  for  the  accelerating 
Lorentz  force  solutions.  The  variables  are  sampled  as  documented  above  for  the 
mass-averaging  of  the  regular  and  Mach  reflections,  but  a  different  procedure  is  used 
for  the  MGD  solutions.  The  variables  for  the  MGD  solutions  are  sampled  directly 
behind  the  normal  shock,  well  before  the  solutions  reach  their  maximum  static  pres¬ 
sure.  By  sampling  close  to  the  shock,  the  downstream  characteristics  which  vary 
from  interaction  parameter  to  interaction  parameter  do  not  effect  the  results.  The 
compression  ratio  is  smaller  than  comparable  RR  and  MR  values  because  the  region 
of  high  pressure  is  located  further  downstream  of  the  shock  system.  The  26°  MR  and 
RR  experience  the  highest  static  pressure  directly  behind  the  shock.  The  maximum 
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temperature  can  be  compared  for  all  solutions,  but  the  26°  RR  Tmax  is  larger  than 
theory  would  predict  due  to  the  secondary  MR. 


Table  4.1  Mass- averaged  variable  ratios  through  the  shock,  interactions  (6  =  26°) 


Solution 

Pbs/Pi 

Tbs/Tis 

Tmaxi  K 

Ptbs  / Pti 

Ttbs/Tti 

RR 

22.18 

3.17 

437.13 

0.283 

1.0 

MR 

19.68 

3.63 

427.19 

0.249 

1.0 

Q  =  1.0 

13.81 

3.58 

1114.94 

0.346 

1.21 

Q  =  3.0 

12.78 

3.74 

4981.32 

0.491 

1.47 

Q  =  5.0 

11.35 

3.68 

1852.67 

0.614 

1.61 

Q  =  7.5 

10.64 

3.60 

1444.27 

0.743 

1.72 

Q  =  10.0 

9.47 

3.48 

1804.07 

0.838 

1.78 

In  Table  4.1,  as  the  magnitude  of  the  magnetic  field  increases,  the  compression 
ratio  decreases.  This  parameter  is  highly  affected  by  the  sampling  location  and  this 
trend  is  only  true  of  the  compression  ratios  directly  behind  the  shock.  Table  4.2 
presents  a  comparison  of  the  pressure  and  temperature  ratios  from  the  inlet  to  the 
point  of  maximum  static  pressure  for  the  same  solutions  presented  in  Table  4.1. 
The  percentage  the  MR  moves  as  a  function  of  the  initial  distance  between  the  26° 
Mach  and  regular  reflections  is  also  included  in  the  table.  The  static  pressure  ratios 
through  the  shock  system  presented  in  Table  4.2  show  the  benefit  of  the  accelerating 
Lorentz  force. 


Table  4.2  Pressure  and  temperature  ratios  through  the  shock  interactions  taken 
at  the  location  of  pmax 


Solution 

Pbs/Pi 

Tbs/Tis 

%  Movement  of  the  MR 

RR 

44.4 

4.23 

100.0% 

MR 

28.5 

5.79 

0.0% 

Q=  1.0 

34.47 

6.95 

52.45% 

Q  =  3.0 

41.01 

7.98 

96.54% 

Q  =  5.0 

45.03 

9.89 

109.3% 

Q  =  7.5 

51.49 

8.54 

127.74% 

Q  =  10.0 

53.33 

10.38 

136.96% 
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The  26°  Mach  reflection  has  a  64.19%  lower  maximum  compression  ratio  than 
the  26°  regular  reflection.  Application  of  the  accelerating  Lorentz  force  for  Q  =  1.0 
increases  the  maximum  compression  to  77.64%  of  the  RR  value.  An  increase  in 
the  interaction  parameter  to  10.0  raises  the  maximum  compression  ratio  beyond  the 
value  of  the  RR  by  20.11%.  This  is  a  result  of  the  stronger  shock  system  generated 
by  the  local  M  =  5.44.  This  proves  the  application  of  MGD  is  capable  of  restoring 
and  even  increasing  the  poor  compression  ratio  of  a  MR. 

In  Table  4.1,  the  temperature  ratio  and  corresponding  maximum  temperature 
increase  until  Q  =  3.0  and  then  they  decrease.  When  Q  =  10.0,  the  temperature 
rise  is  less  than  the  Q  —  1.0  temperature  rise.  Table  4.2  does  not  duplicate  this 
trend,  with  the  temperature  ratio  increasing  from  Q  =  1.0  to  Q  —  10.0.  The 
discrepancy  between  the  two  sets  of  data  is  a  byproduct  of  the  sampling  location. 
The  region  where  the  maximum  temperature  occurs  remains  in  roughly  the  same 
location  relative  to  the  shock  as  the  interaction  parameter  increases.  The  location 
of  maximum  static  pressure,  where  Table  4.2’s  data  originates,  moves  downstream 
relative  to  the  shock  as  the  interaction  parameter  increases.  The  data  for  Table  4.2 
still  captures  larger  temperature  ratios  due  to  the  downstream  sampling. 

The  total  pressure  losses  through  the  shock  decrease  with  increases  in  the 
interaction  parameter.  This  trend  is  due  to  both  the  reduction  in  the  size  of  the  MR 
and  the  added  energy  due  to  the  energy  interaction  term  in  Eqn.  2.34.  The  Q  =  10.0 
solution’s  Mach  stem  height  is  5%  of  the  26°  MR  stem  height.  This  reduction  equates 
to  a  smaller  normal  shock  and  larger  oblique  shocks  making  up  the  shock  system.  At 
a  set  Mach  number,  an  oblique  shock  has  a  lower  total  pressure  loss  than  a  normal 
shock  of  the  same  height.  The  additional  increases  in  total  pressure  recovery  due  to 
the  energy  interaction  term  make  it  possible  to  actually  increase  the  total  pressure 
through  the  inlet,  if  enough  energy  is  added  to  the  system. 

The  total  temperature  ratio  in  Table  4.1  measure  how  much  energy  is  added 
to  the  system.  An  inlet  with  no  added  energy  has  a  total  temperature  ratio  of  1.0. 
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When  the  accelerating  Lorentz  force  is  applied,  the  ratio  increases  by  as  much  as 
75%  for  the  Q  =  10.0  MGD  force.  Even  with  this  increase,  the  application  of  an 
accelerating  Lorentz  force  has  increased  the  compression  ratio  across  the  MR  shock 
system  and  the  total  pressure  recovery,  even  increasing  both  values  above  the  26° 
values  when  Q  >  5.0. 

4-4-3  Electric  Field  Magnitude.  Increasing  the  magnitude  of  the  applied 
electric  held  is  an  alternative  approach  to  creating  a  larger  accelerating  Lorentz  force 
(Fig.  3.6)  to  transition  the  MR.  Eq  is  the  non-dimensional  magnitude  of  the  applied 
electric  held.  The  electric  held  is  non-dimensionalized  by  Vo  and  Bq,  so  the  nominal 
/C  =  1.2  equates  to  an  applied  electric  held  of  2032.23  v/m  when  Q  =  1.0.  Increasing 
E0  above  the  nominal  value  proved  to  be  detrimental  by  overheating  the  local  how 
quicker  than  the  how  could  be  accelerated,  lowering  the  local  Mach  number.  The 
induced  heating  of  the  how  is  linear  with  Q,  but  increases  by  the  square  of  E0. 
Figure  4.23  shows  the  converged  solution  when  /C  =  1.8,  an  increase  of  50%.  As 
expected,  the  heat  addition  to  the  shock  system  due  to  the  applied  force  is  too 
great.  Even  with  the  accelerated  howheld,  the  Mach  reflection  moves  upstream. 
This  is  the  beginning  of  inlet  unstart.  This  movement  raises  the  total  pressure  losses 
and  decreases  the  compression  ratio. 

With  /C  =  2.4,  an  applied  electric  held  increase  of  100%,  the  accelerating 
Lorentz  force  causes  the  inlet  to  unstart  as  shown  in  Fig.  4.24.  Therefore,  /C  =  1.2 
was  used  in  all  solutions  presented  above. 

4-5  Magnetogas  dynamic  Grid  Convergence 

Throughout  this  section,  it  was  assumed  that  magnetogasdynamic  effects  are 
independent  of  grid  resolution  and  follow  the  same  trends  as  the  fluid  interactions  in 
Figures  4.3  and  4.6.  The  coarse  mesh  was  used  to  investigate  preliminary  effects  for 
the  many  parameters  in  this  research,  but  the  actual  effects  due  to  the  MGD  force 
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Figure  4.23  Pressure  contours  for  Q  =  1.0,  /C  =  1.8 


Figure  4.24  Pressure  contours  for  Q  =  1.0,  /C  =  2.4 
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were  decreased  when  the  solution  was  resolved  on  the  fine  mesh.  Figure  4.25  shows 
the  converged  solutions  after  an  accelerating  Lorentz  force  is  applied  with  Q  =  7.5 
to  the  26°  MR.  The  MGD  force  on  the  bottom  has  transitioned  the  coarse  mesh  MR 
to  a  RR.  On  top,  the  same  force  has  pushed  back  and  reduced  the  normal  shock  on 
the  fine  mesh,  but  a  MR  remains.  When  the  accelerating  MGD  force  is  removed,  the 
coarse  mesh  transitions  to  a  26°  RR,  whereas  the  fine  mesh  solution  converges  to  the 
26°  MR.  Both  grids  illustrate  the  effect  of  magnetogasdynamics,  but  the  accelerating 
Lorentz  force  does  not  seem  to  be  grid  independent.  A  detailed  magnetogasdynamic 
grid  resolution  study  is  needed  to  determine  grid  dependency. 


Figure  4.25  Accelerating  Lorentz  Force  with  Q  =  7.5,  coarse  mesh  on  bottom  and 
fine  mesh  on  top,  26°  x  26° 
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5.  Conclusions 


The  sensitivity  of  regular  reflections  and  Mach  reflections  in  a  double-fin  inlet  have 
been  examined  numerically  at  =  4.96.  The  solutions  have  been  used  as  the  ini¬ 
tial  conditions  for  numerical  computations  with  accelerating  and  decelerating  Lorentz 
forces.  These  forces  have  been  applied  with  the  intent  of  determining  the  sensitiv¬ 
ity  of  the  Mach  and  regular  reflections  in  the  dual-solution  domain,  to  transition. 
The  parameters  governing  application  of  these  magnetogasdynamic  forces  have  been 
analyzed  to  determine  the  optimal  application. 

The  regular  reflection  has  been  shown  to  be  very  sensitive  to  a  decelerating 
Lorentz  force.  Continuous  application  of  this  force  produced  inlet  unstart.  A  deceler¬ 
ating  force  applied  for  2.36  x  10“5  s,  but  not  longer  than  8.267  x  10~4  s,  transitioned 
the  regular  reflection  to  a  Mach  reflection  after  the  force  had  been  removed.  Contin¬ 
uing  the  application  beyond  this  time  period  drove  the  normal  shock  upstream  and 
eventually  expelled  the  shock  system  from  the  inlet.  The  careful  use  of  MGD  forces 
in  the  vicinity  of  a  RR  is  needed  to  prevent  the  detrimental  effects  of  transition  to 
a  MR  or  inlet  unstart. 

The  accelerating  Lorentz  force  has  been  shown  to  be  very  sensitive  to  the  speci¬ 
fication  of  the  conductivity.  The  location  of  the  force  determined  by  the  conductivity 
affects  the  location  of  the  Mach  reflection  and  the  additional  effects  to  the  flowfield. 
The  accelerating  magnetogasdynamic  force  as  presently  applied  minimizes,  but  does 
not  eliminate,  the  Mach  reflection.  This  yields  increased  total  pressure  recovery 
by  the  shock  system,  increased  static  pressure  compression,  and  downstream  move¬ 
ment  of  Mach  stem  to  a  more  stable  position.  Depending  on  the  magnitude  of  the 
magnetic  field,  the  accelerating  Lorentz  force  proved  capable  of  increasing  the  total 
pressure  recovery  and  the  static  pressure  compression  beyond  the  regular  reflection 
values.  This  application  resulted  in  a  minimal  Mach  reflection  further  downstream 
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than  the  regular  reflection.  Therefore,  the  accelerating  force  has  shown  the  capability 
to  prevent  inlet  unstart  and  increase  inlet  efficiency. 

Future  research  on  the  use  of  MGD  to  control  Mach  reflections  is  needed  in 
three  areas.  A  magnetogasdynamic  grid  resolution  study  must  be  performed  to  deter¬ 
mine  grid  dependency.  The  MGD  forces  have  been  shown  to  remain  grid  dependent 
for  mesh  resolutions  beyond  which  fluid  interactions  have  been  shown  to  be  grid  in¬ 
dependent.  This  research  will  include  an  examination  of  basic  magnetogasdynamic 
relations  on  grids  of  varying  resolution.  The  viscous  effects  of  laminar  and  turbulent 
boundary  layers  on  the  flowfield  with  accelerating  and  decelerating  Lorentz  forces 
needs  to  be  investigated.  The  added  viscosity  could  effect  the  results  of  the  applied 
force.  In  addition,  the  MGD  code  needs  to  be  modified  to  calculate  the  ratio  of 
specific  heats  based  on  the  temperature.  This  modification  could  be  as  complex  as  a 
multi-species  reacting  flow  calculation,  or  as  simple  as  a  table  look-up  based  on  the 
temperature. 
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